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Introduction 


The  past  five  years  of  this  project  have  been  devoted  to  development 
of  qualitative  theories  of  the  fundamental  electronic  structure  of  pure  and 
impure  semiconductors.  The  basic  theoretical  tools  and  also  conclusions 
relating  to  chiefly  photoemission  type  experiments  and  to  some  extent 
spectroscopy  of  pure  II-VI  and  III-V  compounds  are  presented  in  Chapters  1 
and  2  of  this  report .  Basic  complete  work  relating  to  adsorbates  and  their 
spectra  are  reported  in  Chapter  3.  Several  related  studies  are  being  com¬ 
pleted  at  this  point  and  several  articles  relating  to  these  studies  are 
being  prepared  for  publication. 


\ec\'/ 


Ckt^pter  1,  _ .  ■ 

y  ^  '■t  "J-f  ~^/  c  ,v.^6'c;  1.;^^  O 

Introduction 


This  investlgfjtion  has  the  purpose  of  performing 
self-consistent  energy  band  calculations  on  some  of  the  II-VI 
compounds,  such  as  cadmium  sulfide,  zinc  oxide,  and  zinc  sulfide. 

There  are  many  practical  end  thecretical  reasons 
for  interest  in  these  materials. 


ra.f  CdS). 
C'U  3  or  C 


These  ra terials  have  applications  as  phosphors  (ZnS 
as  infrj-'red  detecrors  (ZnS),  in  photovoltaic  cells  (Cd.l- 
il'-ldVe  heherojunctior.n )  ,  in  catteries  (ZnO),  as  FETs 


2 


(Cc'S),  in  heteroj 
amplifiers  (CdS), 
■j’.hey 

v:ell,  as  e'/amolcs 
the  covalent  III- 
1  t;S i;  la.  tors  .  They 
used  on  th;e  I'ght 


u.iction  lasers  (CdS),  and  even  as  acoustic 

arc  of  considerable  tiieoretical  Interest  as 
of  v.'lde  band  gap  .solids,  intermediate  between 
V  serr.iccnductors  and  the  more  io’.lc  I-VII 
serve  as  a  test  of  calculationr.l  methods  previously 
or  III-V  compounds  and  ether  compoundG  si-Tiplsr 


V]  cr;  :rou.idK , 


then  these 


In  particular,  the  occupied  d~levels,  ivhlch  lie  in  energy  near  the 
top  of  the  valence  levels,  complicate  the  picture.  The  positions  of 
the  d-bands  have  not  been  correctly  predicted  in  previous  ab  Initio 
calculations,  and  in  fact  some  methods  fail  qualitatively  by  not 
obtaining  the  proper  ordering  of  the  valence  s,p  and  d  levels. 

The  II-VI  compounds  have  also  been  of  theoretical  interest 
as  excitonlc  systems.  High  densities  of  such  electron-bole  complexes 
have  been  produced  and  observed  in  CdS .  The  electron-hole  liquid  has 
been  observed  and  studied  in  CdS  as  well. 

In  addition,  some  very  interesting  rragnetic  properties 
in  CuCl  and  CdS  ha.ve  recently  been  observed.  This  is  the  background 
of  these  experimental  studies.  Anomalous  diamagnetism  has  been 
observed  in  CuCl  under  rapid  pressure  and  temperature  cycling.  In 
some  samples,  up  to  50/^  of  the  magnetic  flux  is  excluded.  This  occurs 
at  high  temperatures  (I50  k!  )  and  cannot  be  explained  by  ordinary 
diamagnetism.  The  only  currently  known  phenomenon  that  can  explain 
this  is  superconductivity. 


9 


It'S  clear  that  this  v;ould  have  to  be  superconductivity 
of  different  origin.  Acoustic  phonons,  the  known  mechanism, 
are  too  low  in  energy  to  bind  electrons  at  such  high  temperatures. 
Models  using  some  other  form  of  Interaction  to  pair  electrons 
seem  to  be  necessary.  Some  form  of  excltonic  binding  is  proposed 
in  many  of  these  models.  Cne  of  the  most  well-known  such  models, 
the  Allender-Eray-pardeen  model,  requires  a  metal  and  a 
semiconductor  to  exist  in  very  close  proximity.  Since  CuCl  can 


disproportionate  via  the  reaction  2CuCl 


Cu  +  CuCl  ,  this  model 
2 


has  r.nncffleJ:  to  r.'triy.  Howver,  cudr.i\:.r.  sulfide,  which  has  been 
observed  to  have  only  one  valence  state  and  is  not  believed  to 
disproportionate,  has  shown  similar  anomalous  diamagnetism.  The 
lav;  of  Tvarslmony  is,  therefore,  against  the  metal-semiconductor 
sandv.’ich  idea.  Also,  in  both  CuCl  and  CdS,  impurities  clearly 
play  a  major  role.  Sufficiently  clean  samples  display  no 
interesting  behaviour  at  all.  And  to  make  life  even  iriore 
Interesting,  cadmium  sulfide  also  becomes  a  ferromagnet  under  e 


high  (^1-0  kllognuss)  applied  field. 


{ 

Another  iccdel,  proposed  by  Bishop  and  Overhauser  (  ),  involves 

the  interaction  mediated  by  the  optical  phonons.  This  mechanism 
has  deep  attractive  potential  vjells,  about  10  mev  in  depth  for 
a  spacing  of  250  angstroms.  This  suggests  the  possibility  of 
bound  states,  and  the  pairing  is  considerably  stronger,  around 
a  factor  of  ten,  than  the  pairing  from  acoustic  phonons  in 
the  typical  Cooper  pair  in  the  highest  temperature  superconductors 
I'no'.rn,  at  app-rorimately  20  degrees  Kelvin,  This  idea  is  even  more 
intriguing  in  the  light  of  the  fact  that  both  CdS  and  CuCl  are 
strongly  polnr  compounds.  This  suggest  that  these  phenomena 
v;Guld  not  be  seen  in  the  less  polar  III-V  compounds  and  in 
elemental  semiconductors  suc’n  as  Ge  and  Si.  They’ve  been 
ir.vestlgated  most  thoroughly  and  indeed,  such  effects  have  not 


beer,  observed. 


Chapter  Cne 


The  calculations  discussed  here  are  based  on 

Hartree-Fock  theory,  the  most  common  approximation  to 

the  exact  non-relativistic  theory  of  a  many-eleotron  system, 

1*11  discuss  the  theory  and  its  applications  to  crystalline 

solids.  The  specialization  of  Hartree-Fock  theory  used 

in  these  calculations,  the  method  of  local  orbitals, 

deserves  and  pets  a  chanter  of  its  ovm. 

The  original  prohlei:.  is  finding  exact  solutions 

to  the  Schrodingor  equlation  for  a  nany-electron  system. 

This  is  impossible,  if  the  universe  is  ns  small  as  we  think 

it  is.  Consider  the  wnvefunction ,  vfhloh  is  a  3N  dimensional 

function  if  there  are  N  electron.s.  Dividing  each  axis  into, 

say,  100  units,  in  order  to  numerically  integrate  and 

6n 

differentiate,  we  need  to  record  10  entries.  The  same 
problem  arises  if  we  try  to  describe  the  wavefunction  by 


sets  of  orthogonal  functions. 


6n 

10  is  Just  too  big  a  number. 

Quantum  mechanics  is  agalanst  us.  The  upper  limit  for  the 

2 

rate  of  information  flov;  in  a  computer  of  mass  M  is  Kc  /h 
(h  s=  Planck's  constant,  c  =  speed  of  light),  in  bits  per 
second.  It's  not  hard  to  see  that  macroscopic  physical 
systems  cannot  b.s  exactly  simulated  by  digital  computers. 

I  emphasize  the  v;ord  digital,  since  the  electron  is  an 
excellent  analog  of  itself. 

Such  a  digital  simulation  xvould  be  useless  even 
if  it  could  be  done.  As  suggested,  the  real  solid  simulates 
Itself  perfectly.  The  goal  of  physics  is  to  explain  some 
set  of  occurrences  in  terns  of  a  fe;;,  relatively  simple 
ideas,  rather  than  making  a  model  isomorphic  with  the  universe 
and  just  as  confusing. 

* 

We  begin  to  clear  away  some  complications  by 

using  the  Bcrn-Cppenheimer  approximation.  This  neglects 

any  relationship  betvfeen  the  motion  of  the  electrons  and 

the  motion  of  trie  nuclei.  Inis  a  valid  simplification,  since 

3  5 

the  nuclei  range  from  10  to  10  times  heavier  than  an  electron. 


V/e  can  express  the  K-clectron  v.pvefuncti on 


as  a  function  of  the  electronic  coordinates  alone* 


-  r(,,y  ...  >>  ) 

/y  A  A 


(1.1  ) 


Here  contains  the  position  and  spin  coordinates  of  the 


ith  electron.  We  v;rlte  the  Schrodinger  equation  for  this 


v^avefunction : 

Here,  K  is  the  non-relativistic  p{arniltonian  expressed  in 


atonic  units.  In  this  syste.m  of  units,  =1,  c  =  2,  end 


the  r;as£  of  the  electron  is  or.e-hnir.  The  unit  of  distance 


is  one  bohr  (.52$  angstroms)  and  the  unit  of  energy  is  the 


rydberg  (13.6  ev). 

H  =  H  f  H 

e  r\  tnT 


(1.3) 


//,  =  -  J  y  -f- 1 

-  J.  A  i?  nhf, 

£t  ’  ' 


(1.4) 


j  ^  X  //?  -  i-i  ; 

X 

:i  is  inri 0 pendent  of  the  electron  coordinates,  and  has  been 
n 

assi:!.;f:d  v he  o  constant,  as  per  the  horii-Oppenhe  i  irier 


re  no::  1 t  ion  , 


The  upper  case  characters  refer  to  nuclear 
properties:  Z  is  the  atcaic  nuinber,  H  is  the  nuclear  position 
(  usually  assumed  to  be  fixed,  like  Pimlico)  and  M  is  the 
nuclear  mass.  The  lower  case  characters  describe  properties 
of  the  electrons.  >f < is  the  coordinate  of  the  ith  electron. 

A 

The  atonic  numbers  of  the  atoms  under  consideration 
ere  relatively  low,  and  since  in  any  event  it  is  the  valence 
electrons  that  are  of  primary  interest,  the  non-relativistic 
nature  of  the  Hamiltonian  r.ny  be  acceptable. 

The  fundamental  arproxir.a ticu; ,  r.ov;  applied,  is  the 
independent  particle  model.  In  truth,  in  the  real  solution, 
the  variables  are  not  separable.  V/c  raust  assume  that  they  are, 
or  more  exactly  that  the  true  wavefunotion  can  be  v?ell- 
approxlmated  by  this  roodel.  The  independent  particle  model  says 
that  the  electron  is  acted  on  by  the  average  of  the  other 
electrons.  The  X-electron  v.’avefunction  is  expressed  as  a  product 
of  cne-electron  v.-avefunctions ,  or  as  a  linear  combination  of 
such  products.  A  simple  example  is  the  Hartree  v.’avefxinction : 


(1.5) 


Since  clectronr.  are  ferrier.s :  this  r.-e?^ectron  v.’ovef unction  jr.ust 
be  nntisyratr.etric .  Any  interchange  of  tv;o  f-articles  ir.iist  reverse 
the  sign  of  the  v.-sve function.  The  simplest  v;ay  to  antisymmetriT^e 


is  to  use  the  determinant  of  o  matrix  v.'hose  elements  are 


splnorbitals. 


I  =  (A/If Dot 


DUD  .  ' 


(1.6) 


0(^,)  *  '  ' 


B(^  ) 

A/ 


!:ot  all  such  S^.  are  spin  c irensta tes ,  but  fortunately  v;e  have 
here  closer- sh.:j  1  fiysto:;:'  "hat  can  bo  exp^esred  as  single- 
determiner,  tal  '.sa.vef  unctlenc . 

For  any  arbitrary  ^  ,  kc  can  find  the  cxoectatlon 

V 

value  of  the  Hsriltonian.  It  is  a  sura  of  one-  and  two-  particle 
3  ntegrals : 

<  /  ///  <fj  -~Ux  /V,/,  >  I- Pjhj>0  .7) 

t  D  ~  ^^  0/ tP  (>')  (1.8) 

r  is  a  pcr.mutnticn  operator,  constructed  so  that 


r  /  a  / 


,  The  Ic'-.rs  i:i.i  arc  not  n  nroblem  since  the 


scl  r--C!icr.  y  rr.r.colr:  the  '  r.c If-cxchnnp e  '  . 


The  expectation  value  of  thfi  ilair.lltonlan,^  ^ \  can  be 

N  71/ 

minimized,  by  varying  the  spinorbltals.  It  has  been  shovm  that 

i  !  %,'>  2.  t  ('-9) 

<  t  rt) 


where  E  is  the  energy  of  the  exact  ground  state. 

We  vary  the  orbitals,  but  take  care  to  keep  them 
orthogonal,  for  convenience  and  v;ithout  loss  of  generality.  This 
is  done  by  using  Lagrange  multipliers.  The  objective  is  to  obtain 
the  lov.’est  upper  bound  to  the  exact  non-relativistic  energy. 

We  define  the  functional 


and  require  that  t;ie  variation  of  L  =  0  for  1  =  1,2,...K  . 


0) 


This  yields  a  set  of  i:  coupled  nonlinear 


in tegrodif f erent ial  equations  known  as  the  Kartree-F,ock  equations. 

i  )]  ^  if 


J  and  K  have  the  familiar  Coulomb  and  exchange  operator  forms. 

&.U:) 


^  C  lh~y'/ 


K-e  w  -.i  yv]  (S(x) 

/  \r-y  ! 

These  equations  n.ny  be  written  In  matrix  fori:'; 


(1.12) 
(1.13)  ■ 


F  0  '  -  s  0^ 


(1.14) 


F  Is  the  I’xN  matrix  of  the  Foc>;  oporritor.  For  closed 


shells,  it  is  Kermltian,  so  S  may  he  dia»:onplized  by  a  unitary 
transformation.  Fortunately,  the  Il-VI  compounds  have  closed 


shells.  We  may  write,  F ~  '> 


(1.15) 


Combinins  this  expression  vfith  equation  1.7  gives  the  Hartree-Fcck 
ground  state  energy, 

F.'  =(1/2)7  f  .  (1.16) 

HF  V  ^  ^ 

The  are  the  sijigle-pnrticle  energies,  almost.  Kooprran's 

A 

theorem  states  t'nat  given  a  variatior.ally  stationary  state 

formed  fro’':  N  simultaneous  spinorbitr.l  (pseudo)eigenfunctions 

of  the  K  particle  Fock  opera toi’  ?  ,  the  st«.t(;s  t^nd  ^  , 

/V'v  AJ-J 

formed  by  altering  the  list  of  occupied  eigenfunctions  of  F  by 
one  entry,  ore  stationary  v;ith  respect  to  further  variation  of 
that  spinorbital  whose  occupancy  has  been  altered.  If  we 
neglect  relaxation  effects,  is  the  entire  change  in  system, 


energy  v.hen  an  electron  is  added  or  removed. 


In  order  to  equate  the  C .  ^  with  .single-particle 


er.erglcs,  v:e  must  also  assume  that  relaxation  effects  are  small  - 


occupied  orbita'lf'  change  only  slightly  when  an 


that  the  othri’ 


There  are  such  changes,  since  the  Pock  operator  is  a  functional 
of  all  of  its  ovm  solutions,  and  changing  one  of  them  must 
change  all.  The  question  is,  how  much?  These  relaxation  effects 
vjlll  be  discussed  in  a  later  section. 


The  I'ethod  of  Local  Orbitals 

The  method  of  local  orbitals  is  a  variational 
technique  in  which  solutions  to  the  Eartree-Fock  equations 
are  sought  for  a  s.r:all  subsystem  of  the  system  of  interest, 

V7lth  the  caveat  that  all  such  solutions  found  must  be  contained 
within  the  occuyiecl  Fock  space  of  the  original  system.  Sets  of 
such  local  orbitfls,  which  are  not  orthogonal,  can  then 
be  assigned  to  ewch  subunit  of  the  larger  system  using  the 
translational  invariance  of  the  Fock-Dirac  density  matrix, 
thus  spanning  the  occupied  Hilbert  sp3.ee  of  that  system. 

A  single  noniterative  rotation  of  the  local  orbitals  vilthin 
this  occupied  space  then  yields  the  exact  self-consistent 
solutions  for  the  original  Fock  operator. 

If  ^  0^  ~  h  is  the  sot  of  local 

orbitals  satisfying  the  Eartree-Fock  equations  for  .an  n-clectror. 
system.,  then  f~  ~  ^  ^  (3.1) 

The  Focl;-Dlrac  density  orcratoi'  for  the  equation  is 


/ 


•'J 

p  -  i'.,  ^ 

S  is  the  Inverse  of  the  overlap  matrix  S. 


(3.2) 


J,  =  J  />rJ  ^  M  43.3) 

^  -V 

^  z  t 

Rho  iy)  is  a  projection  operator  ^ 

chosen  such  thatj7^=  for  any  orbital  ^'in  the  occupied  space 

andjy^^O  for  any  orbital  X  the  virti;ial  space.  This  means 

that  for  any  one-electron  operator  L,  the  projection  of  L 


onto  the  Hartree-Fock  manifold  will  satisfy 


(3.^) 


J  X;  ‘3.4) 

for  any  occupied  orbital  ,  and^,  for  any  virtual  orbital/. 


Another  approach,  rigorously  developed  by  Gilbert  (  ), 

tegins  by  introducing  modified  Hartree-Fock  equations  of  the  form 

(  F  ^  (3.5) 

Above,  A  is  an  as  yet  unspecified  Hermltian  operator  and  rho 
is  the  'diagonal'  (  x  equals  x  prirac)  form  of  the  Pock-Dlrac 
density  operator  (  ).  We  see  th^t  rho  is  idempotent,  Kermitian, 

and  projects  onto  occupied  Fock  smce.  liecause  of  this  last 
rrcrerty,  it  can  be  shov;n  that  the  occupied  elg;enfunctionE  of 


/(> 

tho  r.ocUfled  n«J.  t:-oo*-?ccV.  equation  ?clc  er.t3rc;?-y  v.'itr;in  the 
occupied  Fock  space  of  the  unmodified  Hartree-Fock  equations. 

/  _  r  ^  cct.~js.c-J 

So*  y  A  "  I  0  ci^t.rss^'.ie  (3.6) 

There  is  no  Kooprr.an’s  theorem  for  this  modified 
equation,  and  so  the  s  cannot  be  interpreted  as  being  nearly 
the  particle  energies,  as  they  can  in  the  regular  Hartree-Fock 
eq’iations.  The  tota],  energy  of  the  subsystem  still  has  significance 
since  u  iinitary  i.;;  ar^sf  orna  t.l  on  v.'ithin  occupied  Fock  space 
lrav‘:s  the  trace  of  ;'er::ii  t.i operntors  such  as  the  ilamlltonian 

A  ,  .  .  f  A. 

i  .  •  *  .  .a  . 

The  jr-orertics  of  the  /^.a’.;J.ll  vary  uithi  the  choice 
of  A.  ’.'’c  arc  loo'.i?  •  for  furctlnns  that  are  locuilizcd  on  specific 
lattice  sites  of  -a  ci'ystal.  in  order  to  decide  v:hich  operator 
best  fulfills  such  f;  purpose,  wc  need  to  consider  the 
r.ir.irlzation  of  <:  f uncti cjn'ci  >1) I  ^  ^  ^  / 

suijT'Ct  to  the  VO  7' j  a  ti  vona  1  r.cnstra  Ijits  <  =0; 

^U-)  -jc)  -0:  >  r-.o  (3.7) 

■  can  use  i':.--  :'.ei>;ol  c'f  j/irrangc  mult.iplicrs. .  The  integral 


constraints  of  equation  3.7  taken  into  account  by  the 


parameters  ^  ,  and  the  constraints  that  involve  exclicit 

functions  of  X  in  equation  3*7  are  taken  care  of  by  I^grange 


functions  (x)  and 


Vie  obtain 


({  Miej  1  3 


.2,  \lx">  (<xjj}  --<  .i  ))  ~2  (j>lj>~}A?)  (3.8) 

•V-  JM  ^ 

“  / *  ^  1  2  j 'ly  ^  \  cj)  9  h  (x)  ~  o 

H.  J 


Kultiplylng  on  the  left  by  (1-p),  v;e  see 


U-^)y>-o 


j=l ,2, .  .  .M 


(-J’)ijp-  }  f  -0  for  all  7l^ 

end  so  J  ^  (x)  p  (  ~/i  (x) 

Substituting  equation  3.1^  into  3.8,  vie  get 


(3.9) 


(3.30) 


(3.11) 


j’i 


"  I  ?  ij-> 


^  c  S--J 

v.'htch  rnay  be  rfcv:rltten  as  0  fk'^  =  V  ^  !  -y  } 

J 


wr.ere 


(3.12) 


(3.13) 


(3.14) 


/  unitary  transformation  gives 

C  //^y  =  €'!«'> 

rs. 


(3.15) 


Since  v:e  nov:  kno.;  how  to  rnlnirize  any  chosen  Hermltlan 


funstional  ,  v:e  can  now  got  l-acV  to  the  true  problem.  Let 

A  C  -  F,  so  that  £  (  f  y  j  ~  ^  (3J-6) 

Vi’e  can  nov;  ininir;i7.e>^i<Jand  solve  eqn^ition  3*5  at  the  same  time, 

and  so  may  pichTi^i't^Jto  satisfy  physical  reasoning. 

VJe  choose  A  so  that  the  solutions  are  localized  on  a  site, 

usually  by  pichir.r  A  to  be  some  potential  well  centered  on  the 

site.  The  Fock  operator  car.  be  separated  into  F  ,  the  atomic 

s 

Fock  operator  for  the  site  s,  and  U  ,  an  external  potential 

s 

operator  for  the  site  s.  V.e  usually  pick  the  localizing  potential, 

A,  to  be  -U  .  T)‘.c-  local  orbitals  equation  is  nov; 
s 

I  F*  iy,  -  pl  0  .  =  t  6^  (3-17) 

-5  J  J  ^  J  JTx 

This  isn't  an  obvious  improvement  over  the  original 

Hartree-Fock  equation,  but  it  allovis  a  systematic  approximation 

that  simplifies  areatly  the  Hartree-Fcck  problem  without  inducing 

undue  error.  Thr:  left  side  of  equation  3*1?  cab  be  analyzed  by 

the  order  of  the  interslte  overlap.  The  Fock~Dirac  density 

cnerator  is  p  --  Z  T  0  V  0 
y  si  7j  -5  / 

with  V  the  Intrn'n  r.omic  overlap  matrix. 


(3.18) 


V .  ^  =  r  6^^.(yJ  0,  U)  r/ j( 

As  discussed  In  Lovidln  (  ),y  iiiay  be  expanded  in  a  pov;er 


(3.19) 


-/ 


series  in  y  ;  \/  =  -  V  +  . . . 

.  -  ,  «  ■  _  s  *  V 

i,  /  "  '1 


(3.20) 


Expanding  the  local  orbitals  equation  and  ditcarding  second 


order  and  greater  in  intersite  overlap  {  )  gives 

IF,  •*- 1//  -y.  ny,  J 


(3.21) 


and 


^21  '} 
T  II-  -y'l 


vhere  p  =  ^ 

‘r  ""  V 

I  {-^  K  -i— 

^  Tts  ^ 

As  self-consistency  is  approached,  the  last  tvio  terns  of  equation 
3.21  tend  to  cancel.  Therefore,  the  FocJc  space  rotation  changes 


Eloch  v.’avefunctions  into  atcnic-like  v.*svefunctions  that  are 


eigenfunctions  of  the  ator.ic  Fock  operator. 


We  procede  to  localize  a  physically  realistic 


number  of  electrons  on  each  site.  In  cadmium  sulfide,  v?e  choose 
to  have  46  electrons  on  the  cadmium  site  and  IB  on  the  sulfur, 
epproxliratlng  an  ionic  charge  distribution  (  overlap  and  diffuse 
valence  orbitals  picrmlt  cc.mplete  covalency  if  needed,  hovrever). 


V.'e  have  achieved  several  aci! vantages.  The  problem  is 


c  or:  pu ta  t  i  or.a  1 1  y  s  i  r; pi  1  f  i  e  d 


since  the  equation  being  iterated 


does  not  have  to  jortray  the  full  cont'lexi ties  of  the  solid, 
and  since  the  local  orbitals  are  recognizably  modified  atomic 
orbitals  and  aid  in  picturing  chemical  bonding.  This  approach 
has  use  as  a  beginning  in  looKing  at  ancrphous  substances,  which 
have  only  a  short-range  order. 

The  local  orbitals  are  calculated  using  the  program 
called  LOrf.S  vrrltten  by  hunz.  i\  basis  function  expansion  of 
Slater  orbitals  (Sics)  is  used  for  the  radial  part  of  the  local 


orbital  while  a  spherical  harmonic  describes  the  angular  behavior 

B,  /  (r)  -  O.d)  I 

V.’e  usually  take  the  /!/  and  ih  .  from  ragus  et  al  (  )  v»here 


optimized  basis  sets  have  been  computed  for  the  atomic  or  ionic 


system  aprropriace  to  the  solid  in  question.  In  many  cases,  ue 
v.-lll  change  this  basis  set  to  cchicve  localization,  considering 


at  the  same  time  the  total  energy  of  t}ie  subsystem. 


The  expansion  coefficient  C  '•  is  solved  by  the  matrix 


method  of  hoothaan  (  ) . 


v.'e  have  only  used  the  long-range  part  of  the  external  potential 

as  the  localizing  potential.  So,  for  ionic  substances,  the 

electrostatic  potential  v;ill  be  screened  vihile  the  short-range 

effects  of  the  inner  shells  viill  change  the  local  orbitals.  Some 

anions  vjill  not  localize  in  this  situation  since  it  is  the 

Kadelung  part  of  the  potential  that  stabilizes  them.  Since  A  is 

essentially  arbitrary,  we  are  free  to  alter  it  in  vrhatever  v;ay 

.  / 

that  achieves  optirr^l  localization,  and  we  often  use- — ^  jy 
in  equation  3*21  Instead  of 


CHA  FTEli 


o 


Energy  Band  Theory 


The  ?Iartree-Focl:  energy  bands  are  obtained  by  plotting 
the  one-eleotron  energies  derived  in  cq\i.ation  l,i6  against  the 
crystal  ino^rentum  k.  Our  calculations  make  use  of  t’ne  symmetries 
of  the  crystal.  In  the  case  of  these  calculations,  it  is  the 


zinc-blende  crystal  structure. 


Bloch’s  theorem  requires  that  the  v.’nvefunctions  of  the 


.solid  .satisfy 


vj  n  Cl'  c 


/?)  - 


yU  '  t'  J 
hT 


(2.1) 


for  any  lattice  trar.slatjon  vector  R.  The  ci'vst''!!  I'lO.’nentuin  is 

again  k  and  n  is  the  band  index.  V;avefunctlons  that  satisfy 

this  conditicn  fire  Bloch  'wavefunc tlons .  ” 

If  v;c  lock  at  the  set  of  R  ,  we  can  see  that  there 

n 

is  fi.n  clcnen.t  of  R  that  has  a  minl.vum  length,  corresponding  to 

n 

the  fundai-- rtal  lattice  spacing.  Conrsidering  this,  and  the  fact 
ti;nt  :i  is  a  group  under  tr/msla  tior ,  we  can  see  that  there  is 


a  ravji;  k  wt'.en  wt  tra rn'; forrr  into  r: omencum  sraco. 


for'^ ,  there  is  finite  re-vion  of  >io)i-equl  wf^ilent  h-points  - 
the  first  Erillouin  ?-one.This  r:«2/  be  reduced  to  on  '  irreducible 
v/edge  '  by  sy.rnrnetry  considers  t ions . 

The  solution  to  the  Hortree-Fock  equations  is  lengthy 
for  crystals,  and  the  successive  approximation  nature  of  the 
calculations  requires  that  the  equation  be  solved  several  tiroes. 
The  Hartree-Fock  equation  can  be  rewritten  in  terms  of  the. 


Fock-Dirac  density  r.atri>:(  ): 


l’  cf 

vr,.;vc  p  '  '' J'’^  ^  ~ 


(2.2) 


(2.3) 


Jf  a  self-consistent  density  matrix  is  available,  then 
eq-.:atiC/r.  2.2  need  only  be  solved  once  to  obtain  the  eigenvalues 
end  eigenfunctions .  Adams  {  )  and  Gilbert(  )  have  shown  that 

the  density  matrix  (x,x’)  Is  the  same  for  cqirations  2.2  and  >.21 


1  ^-lld  )/'  0^  (2.^0 

J  ■'  ».  r*’“  ^  '  1  r  ^  y-  *  ^ 

So  v;c  have  the  required  density  matrix,  and  can  obtain 


cry.st.a]  wavef unctions  from  the  already  obtained  local  orbitals 
by  CO  structlng  the.  Fock  opera  tor  and  solving  equation  3.21  just 


nr.ee.  The  ('at-nix  oler.entj:  o^  m  are  caiculntcd  to  first  order  in 


ir terp. to.Mic  overlap  ,  consistent  v-.-iti;  the  calcnla t i on  of  the 
local  orbitals. 

The  Bloch  functions  ^„{r)  are  expanded  In  a  basis  set 
in  which  the  basis  functions  teve  the  form 


This  linear  combination  of  basis  functions  technique  is  similar 
to  the  famous  linear  combination  of  atomic  orbitals  technique 
(LCiC)  except  for  the  fact  that  the  are  not  free-atom  orbitals, 
but  are  the  local  orbitals  obtained  previously.  Kote,  the  local 
orbitals  ?ire  occasionally  er.riciied  for  the  case  of  virtual  states. 
The  LTTAS  pro;ira;r.  vrill  not  indicate  basis  fxinctions  a pproj-)ria te 
for  a  virtual  state  if  that  virtual  state  is  of  a  different 
anf'ular  niorr.er.tum  type  than  the  occispied  states;  for  instance, 
e  case  in  v.’hich  the  occupied  states  are  all  s  and  p-like,  while 
the  first  virtual  state  is  d-like.  In  this  case,  basis  states' 
for  the  virtuals  are  added,  consisting  of  spherical  harmonics 
multiplied  by  single  STO's  that  are  chosen  to  have  small  overlap. 

.Since  the  crystal  momenturr.  k  is  a  good  q\:antum  number, 
the  ;c-c;k  s;,!r.e  is  d  i  ngona  1  i/.od  into  cetarnte  S{>aco.s  for  each 
rec j r roca 1  lattice  vector. 


Ihc  inte^i'aj.s  ncce.ssai'y  to  perform 


2^ 

the  j.culc  tion  oro  ii-lndetx-'  h  r  t  ^^^;d  need  only  be  porfor:;;ed 
once.  I'ulticentor  integration  In  done  by  the  Lov.din  c<~f unction 
expansion  method  (  ).  One  of  the  .sites  is  chosen  as  center  and 

all  functions  are  expanded  in  terins  of  spherical  harmonics  centered 
on  that  site.  The  calculations  are  performed  using  the  programs 
KHA'ID  and  KZGXiS  v.'ritten  by  Kunz.  The  ovxtput  consists  of  energies 
and  coefficients  of  the  basis  functions  at  20  selected  points 
of  the  Brlllouin  zone. 


.2  b 

■  CorrcOx tlon  Corrections 

/.s  previously  Kentloned ,  Hartree-Fock  theory  does 
not  give  exact  solutions  to  the  true  many-electron  Hamiltonian. 

Since  we  are  using  a  single  determinental  wavefunction ,  electrons 
are  affected  only  by  the  ?ac;on  field.  Electrons  of  the  same  spin 
have  seme  of  their  true  jair  interactions  taken  into  account, 
tut  there  is  no  pair  corr-ri? tion  at  all  btu’een  two  electrons 
of  opposite  spin.  It  is  obvicus  that  the  Coulomb  force  between 
tvro  electrcr.G  is  i r. depend er.r  of  spin,  and  should  keep  electrons 
from  the  near  neighborhood  o'"  any  other  electron. 

V.'e  say  that  the  motion  of  electrons  of  opposite 
spins  is  uncorrelated,  and  that  electrons  ’with  like  spins  are 
incompletely  correl^ated.  The  correlation  energy  is  \isually, 
defined  as  the  difference  cotv;cen  the  energy  obtained  from 
our  mean-field  J.'artrc<?-Foch  theory  and  the  exact  non-relativlstic 
er.ergy  of  the  system.  Thi.s  exact  enei-gy  could  in  principle  be 
cietermiriod  by  vjslr.g  a  computational  techniq\ie  called  configuration 
in  t  civic  tio.n 


In  this  method,  a  v.avc function  is  used  that  is  a  linear 
combination  of  determinants.  It's  impracticnl  for  all  but  the 
smallest  molecular  and  atonic  systems,  and  is  far  too  time- 
consuralng  to  be  useful  in  solids. 

For  solids  in  general,  and  specifically  in  the  case 
of  the  II-VI  compounds  under  discussion,  vfe  attempt  to  estimate 
the  correlation  correction  to  the  Kartree-Fock  calculation.  The 
correlation  correction  for  the  one-electi-on  eigenvalue  is  the 
error  from  Koopmaru  ’  theorem  (  )  and  from  the  Independent 

pr-.rticle  model . 

According  to  Kooprans *  theorem,  an  orbital  eigenvalue 
is  approximately  the  difference  beyv.'een  the  Hartree-Pock  energy 
of  the  K-electron  system  ivith  the  level  occupied  and  the  Kartree- 
Fcck  energy  of  the  JC-l  electron  system  vjlth  that  same  level 


unoccupied.  These  Hartree-Fock  excitation  energies  are 


6”  --  E  -  E 

At  IJI^ 


(4.1) 


for  states  that  are  occupied  in  the  Ilartrce-Fock  ground  state 


-  r  L 

r'  !=■  H  r- 


(4.2) 


for  virtual  state.'-  of  the  Far  trce-Fock  ground  state. 


v.'herj  v.-t!  to  lt.:i'rc<vo  the  Hartrc-o-.-'ocl-: 


Rp'oroximatlon ,  we  v.ish  to  l\eep  tlie  {general  band  scheme,  so 


(*•)  t  /!(-•) 

i  \,’=  t  t 

V\k. 

t/M*.)  .-U^) 


(4.3) 


(4.4) 


where  the  eigenvalues  are  the  exact  eigenvalues  of  the  system. 

We  assume  that  t'le  correlation  energies  are  a  small  perturbation 


of  the  original  L'artree-Fock  bends,  and  we  can  then  v;rite  E  as 


(L)  (L)  (L) 

E  =  }•;  +  E 


(4.5) 


here,  S  is  the  total  correi-ition  energy  of  a  system  with  L 
c 

electrons.  Equotjov'S  4.3  and  4.^;-  nay  nov;  be  v;rltten  as 

0  (::) 

e  r:  e  +  (K  -  E  )  (4.6) 

nVv  nk  c  c 

0  (N)  (k-l) 

e  r.  e  -}-  (E  -  E  )  (4.?) 

nk  nk  c  c 

A  useful  approxi!  riti on ,  valid  for  nonmetals  v/ith 
valence  band  width  Irss  than  tiio  Hartree-Fock  optical  band  gap, 
was  developed  by  i'antelides  et  al  (  ).  They  shov;ed  that  equations 


4.6  and  4.7  r.ay  be  replaced 

0  (N-1  ) 

F.  =.  E  +  E  (h) 
nk  nk  nk 


(4.8) 


(4.9) 


In  these  equations,  ,  E  ^  (e)  is  the  sslf-cr.ergy  of  an  electron 

nk' 

that  occupies  the  one-electron  Kartree-FocV.  state  hk  in  an 

(K-1  ) 

K-electron  system.  S  (h)  is  the  total  energy  change  in  the 

nk 

remaining  N-1  electrons  when  the  electron  occupying  the  state  in 

(N-1) 

question  is  removed.  E  (h)  is  the  self-energy  of  a  hole. 

nk 

Pantelides  et  al  (  )  have  produced  some  model- 

lrc€;pendent  results  concerning  these  self-energies.  They  find 
that,  the  self-energy  of  holes  is  alvjays  positive,  and  therefore 
the  valence  hands  from  the  original  Hartree-Fock  calculation 
always  move  up  on  the  energy  scale  upon  correlation.  These  self- 
energies  grov:  as  v;e  approach  the  bottom  of  the  valence  bands,  so 
on  calance  the  valence  bands  arc  nari’ov/ed.  The  self-energies  cf 
the  lov;er  electrons  in  the  conduction  band  are  negative,  so  the 
conduction  bands  move  dov;n.  The  models  used  have  only  small 
change  in  the  amount  of  shift  for  the  low-lying  conduction  band, 
end  so  v;e  use  a  rigid  downv.ard  shift  for  these,  V/e  end  up 
shifting  the  valence  bands  higher  and  rigidly  dropping  the 
conduction  b.ands.  This  produces  a  smaller  optical  gap  than  that 
o':.:-. ?ned  by  the  ;;nrtreo-Fock  calculations. 


We  rsorc  nearly  npprc  ';.e{i  ll.t:  exT'Oriniental  gap  by  this  narrov/ing. 

These  self-euergies  v;erc  first  calculated  by 
Toyozawa  (  ).  The  theory  of  this  calculation,  the  electronlc- 

polaron  method,  v;ar.  considerably  further  developed  by  A.  B.  Kunz, 
so  that  it  now  predicts  liole  self~energies  as  well. 

In  the  originol  Kartree-Foch  theory  particles  respond 
only  to  the  average  position  of  the  other  electrons  and  ions. 

This  is  obvicusly  incorrect ,  or  at  least  lnco;:iplete ;  tl:e 
independent  charge  v.'lll  T^oiarizo  its  suri'oundi ngs ,  especially  if 


those  surroundin' G  };avc  sufficient 


tihie  to  re£'pond.  The  electronic- 


polaron  model,  drooso^u  tlic  eoxiov'.c  i'ixr  d  electrons  and  valence 


band  holes  vjlt'n  cuianta 


of  the  polarization  field.  These  quanta  are 


excltons.  In  the  model  we  are  using,  the  excited  states  of  the 

(N) 

crystal  are  simulated  bj"^  a  di spsrsionless  band  of  excltons.  E  (e) 
(N-l ) 

and  E  (h)  are  the  interaction  energies  of  a  bare  electron  and 
nk 

hole  v.'ith  this  field.  These  energies  are  called  polarization 


energies.  This  model  uses  second-order  perturbation  theory  to 


caici^inte  the  self-cncrgi  ns . 


f 


<■ 
»  . 


These  polarization  energies  are  calculated  on  the 
basis  of  a  r.odel  in  which  the  hole  or  electron  is  fixed  in  space. 
This  method,  the  Kott-Littleton  method,  uses  a  perturbative 
approach  to  find  the  induced  dipole  moments  on  all  ions  of  the 
crystal  due  to  the  localized  charges.  This  calculation  takes  into 
account  the  field  from  the  induced  dipoles  -  it  is  self-consistent. 
Since  it  takes  a  finite  time  for  the  crystal  to  respond  to  such  a 
change,  moving  charges  siiould  induce  less  polarization  than  the 
model  static  charge.  Therefore,  this  calculation  should  give  an 
upper  bound  for  the  polarization  energies  of  actual,  mobile 
electrons  ana  holes. 

This  electronic-polaron  model  is  basically  a  long- 
range  scheme.  Only  in  the  limit  of  large  distance  can  vjo  find  these 
changes  by  assuming  that  they  can  be  described  by  dipoles. 

For  short  distances,  the  sliape  of  the  vxavefunctions  plays  a  major 
role,  and  quadrupole  and  higher  multipole  effects  cannot  be 
ignored.  The  clianges  in  the  central  atom  and  its  near  neighbors 
certainly  cannot  be  modeled  by  a  dipole  very  successfully! 


Short-rar..^e  correlation  calculates  the  change  of  the 

nearby  orbitals  v.'hen  an  electron  Is  added  or  removed.  To  find  these 

corrections,  v:e  do  simple  atomic  calculations  of  different 

ionization  states,  supplemented  by  cluster  calculations  in  the  case 

of  negative  ions  vfnere  the  added  electron  extends  over  a  significant 

region  of  space.  In  the  limit  of  zero  overlap  betv;een  atoms,  the 

short  range  correlation  correction  to  the  hole  energy  is 
sr 

E  (h)  =  ^t-Z  (4.10) 

E  C  f 


where  (f  is  the  energy  of  the  level  the  electron  was  removed  from 

and  E  is  the  har'rce-Fock  ionization  energy.  The  ionization 
scf 

energy  is  the  difference  betv:een  the  nartree-Fock  energy  of  the 
system  v/ith  the  level  occupied  and  the  Hartree-P’oek  energy  with 
that  same  level  full.  The  short-range  energies  so  calculated 
are  again  upper  bour.r's,  since  electrons  and  holes  are  not  as 

to 

localized  as  this  model  potrays  them. 

The  Kartree-Fock  calculation  and  the  following 
correlation  corrections  give  a  set  of  one-electron  energy  bands. 
They  must  nov;  be  co'rared  to  experiment  to  test  the  adequacy  of 


our  f!  pprox ^  on  t  ion  , 


CADMIUM  SULFIDE 


c 


At  room  temperature  and  standard  pressure  CdS 
tnllizeo  in  the  zinc-blencic  Lattice  v.’ith  a  lattice  constant 


of  5*8^^  angstroms  (10.99  atomic  units)  and  in  the  hexagonal 
ferr.  v.'ith  lattice  constants  of  4.13^3  and  6.7234  angstroms  ( 
7 . 31 3^^  and  12.705^^  atomic  units).  The  zinc-blcndc  form,  the 


subject  of  these  calculations,  is  composed  of  t\70  interpenetrating 

++ 

face-centered  cubic  sublattices  occupied  by  Cd  and  S  ions, 
disriaced  relative  to  each  other  by  1/4  of  the  diagonal  of  the 


U!'iit  cube.  The  symmetry  properties  of  the  zinc-blende  lattice 
2 

(  sr^ce  grout)  T  )  hove  been  discussed  by  I'anmcnter  (  ),  and  in 

d 


folloviing  discussl( 


ro'io’aort,  Smoluchowst.l  and  Wlgncr  (  )  v;lll  be  used. 


We  first  examine  the  energy  states  at  the  Prillouin 


10 

zone  center,  the  point.  Cadmium  ions  contribute  a  filled  (4d) 

shell  to  the  valence  energy  region,  and  sulfur  ions  contribute 

2  6 

a  (3s)  (3p)  configuration.  The  p  and  d  levels  hybridize  to 

some  extent,  since  they  are  separated  by  less  than  nine  ev;  still, 

the  valence  bands  are  not  strongly  hybridized.  The  top  valence 

band  is  predominantly  p-like  (  ever  75%)  and  the  next  lower  band 

is  predominantly  d-type,  again  over  75%^  We  expect  the  lowest 

condiictjon  band  to  be  derived  from  the  cadmium  5®  levels,  as  is 

usval  in  co;;:pour.:is  that  have  any  ionic  character.  In  the  zinc- 

blende  structure,  the  crystal  field  splits  the  fivefold  degenerate 

d  levels  into  a  tricly  degenerate  state  and  a  doubly  degenerate 

1 

state.  The  threefold  degenerate  p  states  stay  degenerate  under 

12 

this  crystal  field  and  transform  like  .  The  zinc-blende  lattice 

15 

has  no  Inversion  symmetry,  so  the  bands  at  the  gamma  point  need 
not  Ystve  a  definite  parity.  The  conduction  band  contains 

•5 

both  p  and  d  contributions,  in  fact. 

Several  previous  calculations  have  been  cone.  The  valence 
.  ■;al]y  ft;;.  Upper  ones,  find  the  lov;or  conduction  bands 


i  ar 


esT 


hnve  been  v.'o?  T -described  by  varioui;  pseudopotential  calculations. 
That  is,  they  accord  with  current  exp. ..t,  »;i, .  -’•or,  since 

the  pseudopotential  method  is  basically  a  parametrization  scheme, 
in  which  the  parameters  are  determined  by  experiment,  little 
critical  new  Information  can  be  obtained.  If  the  experimental 
evidence  has  been  misunderstood  or  misinterpreted,  the  ; 
pseudopotential  v;ill  simply  predict  the  mistakes  or  misunder¬ 
standings  that  it  sprang  from.  A  pseudopotential  calculation 
is  an  aid  to  undea standing,  but  it  tends  to  not  be  falsifiable. 
Such  calculations,  such  as  those  by  Cohen  and  Eebgstresser  (  ) 

also  depend  or.  the  validity  of  tlie  cancellation  theorem,  v;hlch 
is  not  exact  and  v.’hich  is  moh  less  valuable  for  systems  that 
have  localized  states  that  ore  closely  comparable  in  energy  with 
the  valence  states  -  here,  the  cadmium  d  states  are  the  ones 
not  easily  described  because  of  their  local  nature. 

A  related  but  more  theoretically  rigorous  technique, 
the  orthogonalized  plane  wave  method,  has  some  of  the  same  errors. 
In  this  method,  the  valence  and  conduction  states  are  described 


by  plane  v.’aves  ort’nogonalized  to  the  core  states,  taken  as  constant. 


in  these  previous  OF,v  calculations ,  ?is  performod  ty  £uv<owa  and 
Stukel  {  ),  Euweir-a  and  Collins  -=c  al  (  )  and  by  Stukel  et  al  (  ) 

there  have  been  errors  of  up  to  three  ev  in  the  p-llke  valence 
and  conduction  I'ands  in  compounds  made  of  first-roiv  atoms. 

Also,  the  4d  states  are  misplaced  by  about  a  rydberg  -  probably 
related  to  their  core-like  nature.  The  level  ordering  of  the 
cadmium  and  the  sulfur  s  bands  seems  to  be  reversed  from 
that  measured  by  photoemission  by  about  6'ev.  Altogether  it 
is  easy  tc  see  that  there  is  a  need  for  a  first-principles 
calculatiov;  that  correctly  predicts  the  major  features  of 
the  hand  structures  of  these  r.a terials . 

Therefore,  the  first  all-fslectron,  self-consistent  ' 
noner.pirical  band  calculation  usinr  nonlocal  exchange  has  been 
performed  on  CdS.  The  emthods  of  calculation  are  as  discussed 
earlier  in  this  v;ork.  It  should  be  mentioned  that  a  selfrconsistent 


band  calculation  has  been  performed  by  Zunger  and  Freeman  (  ) 

using  H  local  density  approximation  for  the  exchange  potential. 


Their  results,  as  v;ell  as  earlier  ones,  will  be  discussed  in  the 


( 
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d  ly  Tessnan  et  nl  (  )  were  then  used,  for 


Cd  nnd  S  along  with  the  electric.  :  -^ant  of  6.32 

to  calculate  the  polarization  energies  uslt.g  i.;-  Kott-Littleton 

method.  We  obtained  values  of  0.1866  ry  and  0.2534  ry  were 

ob-t-er-ined  as  the  polarization  energies  associated  with  the  S 
++ 

and  Cd  holes  respectively,  peing  less  tightly  bound,  the 
sulfur  anions  polarize  to  a  greater  extent  around  a  hole  at  the? 

-r-»- 

Cd  cation  site  than  do  cadmium,  ions  around  a  sulfur  hole, 

producing  a  larger  polarization  energy  for  the  cation.  We  also 

calculated  s-iort-range  relaxation  energies  for  the  states  of 

sr  — 

interest  in  the  valence  region;  these  are  E  (S  )  =  0,07824  ry, 

3p 

sr  sr  ++ 

E  (S  )  =  0.06302  ry ,  and  E  (Cd  )  0.116  ry .  These  corrections 

3«  4s 

^;ere  added  to  the  lia.miltonian  matrix  v.’hich  was  then  rediagonalized 
to  give  the  correlated  valence  bands.  Conduction  bands  were 
computed  by  shifting  the  Hartrce-Eock  conduction  bands  by  -,2584  ry, 
the  polarize tion  energy  of  a  conduction  band  electron  on  a 

+-r 

Cd  site. 

Shown  in  figure  4.1  arc  the  coi'related  energy  bands 
of  CcuS.  The  calculated  band  sti^oture  shows  cadmium  sulfide  to  be 


Figure  4.1  Correlated  energy  bands  of  CdS  for  the  normal 


lattice  constant  of  5*^81  angstroms  (  10.994  au). 


p  direct  r-rd  ffsp  semiconductor  'with  thei  t’:eD  Pt  the  r;oint. 

This  conclusion  is  in  agreement  \>jlth  the  previous  bend  calculations. 

Three  rr.aior  non-overlepping  regions  constitute  the  valence  bands 

in  this  system.  The  3  3s  derived  band  lies  19  eV  belov;  the 

valence  band  edge  end  is  about  1.4  eV  i^ide.  The  next  region 

++ 

is  primarily  derived  from  the  Cd  4d,  and  lies  12  eV  belov;  the 
top  of  the  valence  band.  T’nis  band  is  about  one  eV  wide.  The 
uppermost  valence  band  is  primarily  derived  frorri  S  3P  levels, 
and  is  about  3.9  eV  v;ide.  The  lowest  conduction  band  is  s-like 


and  is  derived  fro-n  the  Cd  5s, 

3 

3  s  ar)d  S 

4s 

levels . 

The  band  gap  is  found 

to 

be  direct 

and 

eq\i.al  to 

7.1  eV. 

T''e  optical  value  of  the  gap  ,  '£ 

— 

2.55  eV  ( 

). 

cbtn ined 

from 

exr 3r  1  rr.ent ,  is  in  serious  disagreor.ent.  The  correlation  model 
here  used,. the  electronic-rolaron  model  and  it  limit, 

Kotz-I ittleton  theory.  Is  adapted  to  insulators  and  ignores 
shcrt-rar.;:c  polar i/.a t ion  effects.  Such  effects  should  be  small 
in  a  tor  s  and  in  sy.stems  vrhere  the  local  orbitals  are  only  slightly 
pc'^tur  bo  i  I'ror'  tiio  free  a  to.'^  oi-  lor.,  but  in  poliir  semiconductors 


5  •;  3h  ;  s  Pd  ■  ,  thcr 


."ua rar.Pec  thal.  the  local  orbitals  are 


is  Tin 


Figure  4.3  XPS  spectruia  for  CdS  according  to  Ley  et  al  (  ). 

I  and  I  describe  fine  structure  of  the  highest  valence-band 
1  2 

peak.  S  is  a  shoulder  on  the  high-binding  side  of  this  peak. 

1 


II  represents  the  second  valence-band  peak. 


( 

that  slrrilar  to  the  fi*ee  ions.  In  this  class  of  compounds,  there 
is  substantial  charge  density  in  the  Internuclear  region. 

It  has  been  argued  that  such  effects  act  to  r.arrov;  the  band  gap. 

The  density  of  states  for  CdS  as  calculated  from  the 
original  band  structure  calculation  is  given  in  figure  4.2.  For 
coT.rarison,  figure  4.3  shows  the  corrected  XPS  spectrum  obtained 
by  Ley  et  al  (  ).  The  ex  perl  mental  evidence  clearly  shovfs  the 

tv:o  peaks  in  the  density  of  states  of  the  upper  valence  band. 

The  positions  of  these  peaks  are  also  in  reasonable  agreement 
■with  experiment.  The  upper  peak  lies  1.6  cV  belov/  the  top  of 
the  valence  band  by  Ley's  measurements,  -while  our  calculation 
rives  a  peak  at  approximately  1.4  eV.  Experimentally,  the  second 
peak  is  at  4.1  eV,  while  this  calculation  has  a  peak  at  5*4  eV. 

It  also  seems  possible  to  identify  the  shoulder  of  the  upper  peak. 
From  our  band  calculation,  this  shoulder  seems  to  be  at  about 
2.6  eV,  while  Ley"s  measurements  put  it  at  2.1  eV. 

The  position  of  the  d-like  levels  is  also  correctly 
prod ir. ted  in  these  cnlcula tiens .  These  levels,  primarily  formed 
frer.  cf.dmlun  4d  .states,  are  found  to  peak  at9;64  eV  belov; 
the  top  of  t}ie  valci'.ce  band  expci-lmentallly ,  according  to  Ley. 


Figure  4.5  Iniagirury  part  of  the  dielectric  constant  of  CdS 

(  this  calculation) 


Figure  4.6  Experimentally  derived  imsiginary  part  of  dielectric 
constant  in  CdS,  from  Cardona  et  al  (  ). 

Derived  from  ultraviolet  reflection  measurements. 


Cur  C'llcu]'  t  ior.  p  xvos  a  v.-ilue  of  10.-  eV  'nc-lov;  the  tor  of  the 
valence  band,  In  very  Kood  npreement  with  exper i nent .  This  is 
in  sharp  contrast  v;ith  prev'ious  pseudopotential  and  OF'./ 
calculations,  v;hlch  nispiaced  this  band  as  much  as  12  eV  below' 
its  actual  position.  It  should  be  mentioned  that  Ley's 
photoeraission  studies  were  done  on  CdS  in  the  hexasconal  form, 
but  since  the  first  two  shells  around  the  central  atom  ore 
identical  in  the  hc::n£ronal  and  zinc-blende  structure,  the  atom 
is  exposed  to  a  very  similar  potential,  and  no  substantial  change 
in  the  density  of  stfites  Is  expected. 

In  figure  the  joint  dtnu-lt.y  of  states  as  calculated 

from  our  band  structure  is  shown.  In  figure  ^>-.5,  the  imaginary 
part  of  the  dielectric  function  is  calculated  from  the  same 
theoretical  band  structure.  In  figure  4.6,  Cardona  et  al  extract 
the  imaginary  part  of  the  dielectric  constant  from  ultraviolet 
reflection  measurements  on  zlnc-blendc  CdS  (  ).  There  is  deal' 

qualitative  agreement,  especially  jith  the  central  peak  nnd  the 
two  s\).a-r'(;!ys  o  pprox  i  ma  tel  3  and  5  above  that  centred  peak. 

The  tail  hr  lev;  the  central  T'i  f:k  in  energy  Is  exaggerated  due  to 
t  ':u’  f  verer;'.  *.  r  a  te  of  th.c  I'ai  d  gap  in  our  calculations. 


( 

To  sum  up;  these  calculations  correctly  predict  the 
cajor  features  of  the  energy  bands  In  CdS,  ns  determined  from 
optical  and  other  measurements.  Insofar,  this  first  all-electron, 
self-consistent  non-emnirical  method  has  succeeded.  The  d-like 
levels  are  correctly  predicted,  as  vias  not  the  case  v;lth  earlier 
calculations.  Nothing  in  the  pure  bulk  energy  bands  suggests 
an  explanation  of  the  anomalous  diamagnetism,  and  it  would  seem 
that  the  effect  is  not  an  intrinsic  one.  It  may  be  hoped  that 
this  better  understanding  of  the  electronic  structure  of  the 
pure  solid  will  aid  in  the  understanding  of  the  anomalous 
diam^rnetic  state,  perhaps  by  acting  as  necessary  first  step  in 
the  ey.ar.i nation  of  the  properties  of  defects  and  impurities  in 
cadm.iur:  sulfide. 


ZINC  OXIDE 


At  room  temperature  and  standard  pressure  ZnO 
crystallizes  in  the  hexagonal  (  wurtzite)  form  viith  lattice  constan 
of  angstroms  (6.U’30  atomic  units)  and  5»^-93  angstroms  (9.32^ 

atomic  units).  The  zinc-blende  form  with  identjcal  nearest-neighbor 
distance  has  a  lattice  constant  of  ^.595  angstroms  (3.68^4-  atomic 
units)  and  is  the  subject  of  these  calculations.  It  5s  composed  of 
two  interpenetrating  face-centered  cubic  sublattices  occupied  by 
Zn  and  0  ions,  displaced  relative  to  each  other  by  1/4  of  the 
diagonal  of  the  unit  cube, 

V.'e  begin  again  by  examining  the  energy  states  at  the 

10 

Hriilouin  zone  center.  Zinc  ions;  contribute  a  filled  (3d) 

shell  to  the  valence  energy  region,  and  oxygen  ions  contribute  a 
2  6 

(2s)  (2p)  conf igurntion.  In  this  crystal  structure,  the  threefold 
degenerate  r  states  st'-jy  cl e go v, crate  and  transform  lilcc  ’ 


15 


The  crystal  field  splits  the  fivefold  de,jenerato  dlevels  into  p 
trlplv  depcencTp tc  '  state  and  a  doubly  deKcreratc  /'  state. 

'  15  12 

The  zinc  3d  levels  and  the  oxygen  2p  levels  lie  quite  close, 
and  this  calculation  shovis  shov;s  all  eight  valence  bands  lying 
within  a  five  eV  region.  The  zinc-blende  lattice  has  no  inversion 
syr.netry,  so  the  bands  at  the  gamma  point  need  not  have  any 
definite  parity.  This  calculation  shovjs  that  this  upper  valence 
region  exhibits  strong  p-d  hybridization.  The  band  lying 
beneath  this  upper  valence  region  lies  nearly  24  eV  belovr  the 
top  c''  the  valence  band,  and  is  primarily  derived  from  oxygen 
2s  SL:-;tos.. 

In  general,  in  ionic  compounds  such  as  ZnO  we  would  expect 
the  lov;eyt  conduction  band  to  be  primarily  formed  from  the 
zinc  4s  level.  The  actual  calculation  shov;s  that  although  the 
zinc  4s  states  play  a  major  role,  oxygen  2s  and  js  states  actually 
play  ?•  larger  role  in  this  first  conduction  band.  The  next 
cor.d\:ctlon  tarai,  triply  degenerate  at  the  ".one  center,  are  formed 
alr,u>st  entirely  from  zinc  4p  states  .-This  Is  suggested  merely 
by  co'r^^rlng  zinc  with  its  successor  in  the  periodic  table,  gallium. 


Several  previous  calculations  tiave  been  performed  on 
zinc  oxide.  Attempts  using  local  pseudopotential  theory,  such  as 
those  by  Rossler  (  )  and  Bloom  and  Ortenburger  (  )  have  not 

beer,  entirely  successful.  Since  the  pseudoootentinl  theory  depe.rds 
upon  a  cancellation  cf  the  strong  core  part  of  the  potential  by 
the  usual  requirement  that  the  valence  electrons  be  orthogonal 
to  the  core  electrons,  first-rov;  elements  would  seem  to  be 
unsuitable  for  this  approach.  The  point  is  that  valence  p-states 
for  first  row  elements  are  not  required  to  be  orthogonal  to  any 
p-cors  states.  Investigators  have  attempted  to  alleviate  this 
prosier  by  e.rnirical  adjustments  to  the  pseudopotential  s ,  but 
it  hr^sn't  V/orlred  v:ell.  When  nonlocal  pseud.opotf;ntlals  were  used, 
as  in  the  c.alcula  tiers  of  n’nel  iljov/sby  (  ),  a  better  agreement 

with  exrsri.ment  is  achieved.  Even  so,  much  of  the  value  of  the 
pseudopotential  approach  is  lost.  There  is  not  the  same  confidence 
that  the  pse'udopotentia  1  v/ill  retain  its  predltctive  powers  in 
different  compounds,  because  the  physical  rationale  is  weakened. 
The  "odlfied  and  nonlocal  1  dilations  still  suffer  frora  all  the 


ori  rinal  weakness  of  pseudopoter  t  inis ;  reliftnce  on  exyieriment. 
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FUKCTI0X3  FOR  ZIKG  AND  LOCAL  ORBITAL  COEFFICIENTS 
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nccTiGvs  FOB  cxyge:;  and  local  orbital  coefficients 
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Therefore,  vxe  perform  on  all-electron,  nelf-consistent 
noner.pirical  band  colculation  on  ZnO.  The  methods  of  calculation 
are  the  same  as'  used  previously  on  CdS. 

The  calculation  v;as  begun  by  calculating  local  orbitnls 
for  ?:3nc  oxide  v.’ith  the  normal  lattice  constant.  The  zinc  and 
oxygen  basis  sets  of  Bogus  et  al  (  )  were  used  and  vuere  modified 

for  this  calculation.  Two  additional  diffuse  STOs  were  added  to  the 

s  basis  on  zinc.  The  results  of  the  local  orbitals  calculation 
are  listed  in  tables  and  5«2  for  zinc  and  oxygen  ions, 
resroctl vely . 

Hartree-Fock  bands  were  calculated  foi'  20  k  points 

in  the  irreducible  wedge  of  the  Brillouin  zone.  The  usual 

exaggerated  band  gap  appears;  we  proceed  to  apply  the  correlation 

methods  previously  discussed. 

Folarlzabill ties  listed  by  Tessiran  et  al  (  )  were*  then 

used  for  Zn  end  0  along  vjlth  the  optical  dielectric  constant 

cf  •’.'.036  to  calculate  the  polarization  energies  using  the  Kott- 

Littleton  method.  We  obtained  values  of  0.2172  ry  and  0.2364  ry 

— —  ++ 

as  polarization  energies  associated  v?ith  the  0  and  Zn 

holes,  respectively .  We  also  calculated  short-range  relaxation 


igure  5*1  Correlated  energy  bands  of  ZnO  for  the  normal 


lattice  constant  of  ^.595  angstroms  (  8.6S4  au 


energies  for  the  state  of  Interc-.-t  in  the  \'r?;i.onco  roslcn.- 


this 


sr  ++ 

ds  E  (Zn  )  =  0.31^02  ry.  These  corrections  v.ere  added  to  the 

3d 

Kamlltonion  riatrix  vihich  was  then  rediagonalized  to  give  the 


correlated  valence  bands.  Conduction  bands  were  computed  by  shifting 


the  Hartree’Fock  conduction  bands  by  -.2864  ry,  the  polarization 

++ 

energy  of  a  conduction  band  electron  on  a  Zn  site. 


Shown  in  figure  ^.1  are  the  correlated  energy  bands 


of  ZnO.  The  calculated  band  structure  shows  zinc  oxide  to  be 


a  direct  band  gap  semiconductor  v;ith  the  gap  at  the  gamma  point. 


This  conclusion  is  in  agreement  v;ith  previous  band  calculations. 
Two  major  non-overlapping  regions  constitute  the  valence  bands  in 
this  system.  The  highest  region  is  about  4  eV  in  width,  and  is 
divided  into  two  subregions '  viith  almost  no  overlap.  The  higher 
and  wider  of  these  subregions  is  ,  approximately  3  eV  in  width, 

— —  4.^ 

Is  composed  of  O  2p  and  Zn  3d  levels,  vihile  the  Hower  subband 

++ 

is  almost  dispersionless  and  is  of  almost  pure  Zn  3d  character. 

The  lowest  valence  band  is  about  2.3  eV  in  vfidth  and  is  composed 
++ 

of  0  2p  and  Zn  3d  levels. 
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The  band 
the  gap  is  ex 


gar)  is  found  to  be  direct  and  equal  to  10.86  eV 
pcrimcii tally  found  to  be  3.3  eV ,  it  must  be 
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Figure  5*2  Density  of  states  of  ZnO 


( 


Figure  5.3  Corrected  XPS  spectrum  for  ZnO  according 


to  Ley  et  al( 


that  the  conduction  band  structure  here  predicted  is  seriously 
in  error.  Previously  mentioned  potential  problems  in  our  correlation 
model  may  explain  this  discrepancy,  but  it  seems  most  likely  that 
there  is  some  flau  in  the  estimate  of  the  polarization  energies. 

The  density  of  states  for  ZnO  as  calculated  from  our 
band  structure  is  given  in  figure  5»2.  For  comparison,  figure  5*3 
shows  the  corrected  experimental  density  of  states,  derived 
from  x-ray  photoemission  experiments  conducted  by  Ley  et  al  (  ). 

The  d-like  levels  have  been  subtracted  out  from  the  experimental 
data,  since  they  dominate  the  spectra.  The  double  peak  in  the 
density  of  states  of  the  uppermost  band  is  clearly  shovm  in  ■ 
experiment  and  in  our  calculation.  The  d-like  levels  in  our 
/  calculation  have  some  structure,  are  not  smoothed  into  a  single 
.peak  as  in  the  experiment.  The  experimental  peak  is  centered 
8.  Si  eV  belov;  the  top  of  the  valence  band,  according  to  Ley, 

Our  d-co.mplex  is  situated  approximately  5  sV  below  the  top, 
but  the  qualitative  picture,  that  of  a  double-peaked  valence 
band  with  a  very  sharp  d~band  about  3  eV  belovj,  corresponds 


closely  with  exper if^ent . 


Our  calculatiotis  also  {rive  a  band  derived  froiri  0  2s 
states  (not  shovm  on  figure  5*2)  centered  around  24  eV  below 
the  top  of  the  valence  band  -  this  corresponds  to  the  lowest 
peak  in  figure  5*3 

In  addition,  we  have  calculated  the  joint  density  of 
states  (  figure  5*^)  the  imaginary  part  of  the  dielectric 


constant  (  figure  5«5)  fcir  zinc  oxide. 


_£:ure  5..^  Joint  Density  of  States  of  ZnO 
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Figure  5-5"  Imaginary  part  of  the  dielectric  constant  of  ZnO 


CONCLUSIONS 


We've  used  ab  Initio  Hartree-Fock  theory,  as  well  as 
relaxation  and  polarization  correlation  corrections  to  that  theory, 
to  calculate  the  electronic  band  structures  of  cadmium  sulfide  and 
zinc  oxide.  In  the  calculation  of  cadmium  sulfide,  hybridization 
of  the  Cd  ~s,  S  j3  and  S  bs  Is  a  major  factor  in  the  loviest 
conduction  band.  The  calculation  shows  that  CdS  is  a  direct  gap 
semiconductor  with  a  gap  of  7.1  eV  at  the  gamma  point.  This  tends 
to  show  that  the  mechanism  proposed  by  Abrikosov  to  explain  the 
anomalous  dlamagnetisn:  originally  seen  in  CuCl  (  and  discredited 
there  by  Weidman  (  )  )  is  not  applicable  to  CdS,  since  it  requires 

a  small  indirect  gap. 

The  calculated  one-electron  energy  bands  are  compared 
v;ith  the  published  optical  date.  The  valence  bands  are  in  excellent 
agreement  with  photoernission  and  reflection  data,  and,  in  particular 


the 


position  of  the  d-like  band  is  correctly  predicted,  in  contrast 

with  previous  calculations. 

The  calculation  for  zinc  oxide  shows  that  there  is 

very  significant  nixing  in  the  ucper  valence  bands  between  the 
+  + 

C  2p  and  Zn  3d  levels.  It  seens  clear  that  previous  pseudopotential 
calculations  that  could  not  correctly  take  into  account  the  zinc 
d  levels  v:3re  incapable  of  explaing  the  valence  bands.  The  calculation 
gives  a  band  gap  of  10.86  ev,  and  predicts  ZnO  to  be  a  direct  gap 
serai  conductor  with  the  gap  at  the  gamttia  point.  In  the  case  of  both 
CdS  and  ZnC ,  it  should  be  recalled  that  we  are  essentially  just  solving 
Dyscn’s  equation,  and  that  higher  order  diagrams  are  quite  capable 
of  accounting  for  the  difference  between  the  (exaggerated)  gaps 
from  our  calculations  and  the  measured  experimental  gaps.  The  fact 
remains  that  these  calculations  predict  the  valence  structure; 
the  one  electron  bands  for  zinc  oxide  also  accord  well  with  published 
optical  experinents.  Agsln,  the  position  of  the  d-like  bands  is 
correctly  predicted. 

The  anomalous  diairrignstisin  observed  in  CdS  remains 
unexplained.  It  scei.is  to  be  an  exti'lnslc  phenomenon  , 
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The  observation  of  ferromagnetism  upon  an  applied 
field,  as  well  as  the  anomalous  diamagnetism,  is  intriguing. 

It  suggests  that  the  pairing  mechnaism  may  favor  a  state  with  a 
spin  of  one,  as  opposed  to  the  Cooper  pair,  with  a  net  spin  of  zero. 
Each  pair  would  have  an  intrinsic  magnetic  moment.  This  possibility 
has  been  discussed  for  classical  phonon-mediated  superconductivity 
by  ?.  Anderson  and  F.  Eorel{  )  but  seems  not  to  have  been 
observed,  except  possibly  at  extremely  lo\^  temperatures  (  ). 

An  exciton-medlated  pairing  might  well  favor  the  1=1 

or  higher  state.  The  superfluld  state  in  He  is  suggestive. 

3 

AS  yet  the  details  of  the  interactions  are  not  know, 
although  it  does  seem  thnt  chemical  impurities  play  a  role.  It  may 
still  be  possible  to  test  some  of  these  ideas  by  a  phenomenological 
theory  like  the  Ginzburg-Landau  theory.  One  possibility  seems 
Interesting;  in  the  Ginzburg  -Landau  theory,  two  characteristic 
lengths  appear,  the  coherence  length  and  the  penetration  length. 

In  such  a  model  revised  for  a  p-wnve  system,  a  third  length  arises, 
the  characteristic  length  for  the  change  in  spin  direction. 

The  two  lengths  of  Ginzburg-Landau  theory  allow  the 
existence  of  a  distinctive  surface,  v;hlch  may  have  a  lower  energy 


r/\ 
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tr^r.  the  bulk  state,  v;hich  leads  to  the  interesting  and 
technologically  useful  type  II  superconductors.  This  still  exists 
with  three  characteristic  lengths,  but  another  effect  becomes  possible. 

A  second  surface  layer  is  Introduced;  we  may  compare  this  to  the 
earth,  with  crust,  mantle  and  core.  If  the  mantle  is  energetically 
favored,  bubbles  of  a  certain  size  would  be  energetically  favored. 

These  sunerconducting  domains  might  explain  the  very  high  but  finite 
conductivity  seen  in  the  anomalous  states  of  CuCl  and  CdS,  and 
alignment  of  these  domains  might  explain  the  feiromagnetlsra  seen 
in  CdS  under  high  applied  field.  At  the  moment,  this  is  all  speculation. 


We  await  further  experiments. 
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ABSTRACT 

En€;rgy  bands  and  ionisaticr.  energies  for  electrons  in  AlP,  Gap,  and 
GaAs  are  obtained  using  a  nav  netiiod  propo.srd  'ey  Kunz,  ct  al."’  The  valence 
bands  and  ionization  energies  obtained  are  found  to  be  in  good  agreement 
v;ith  experiment.  Tnc  conduction  bands  vary  substantially  from  experiment. 
Finally,  a  method  for  obt£‘,inir.g  better  conduction  bands  is  proposed. 


1 .  IKTrtODUCTIO'J 


A  new  ^  initio  r.ethod  for  obtaining  correlated  enerp.y  bands  vas 
recently  proposed  by  fCunz,  e^t  al..  ^  “Hie  need  for  a  nei;  method  becomes 
apparent  v.-hen  the  inherent  problems  associated  v;ith  current  initio 

methods  are  considered.  Existing  methods  can  be  broken  into  tv;o  classes, 
those  which  use  the  nonlocal  Hartrec-Fock  (H-F)  exchange  potential  and 
those  ivhich  use  some  local  density  approximation  to  the  exchange  or 
exchange-correlaticn  potential.  The  former  class  of  methods  requires 
calculations  which  are  both  lengthy  and  complicated.  Once  such  a  H-F 
calculation  has  beer,  performed,  there  is  still  a  need  to  obtain  correlation 
corrections.  Tneso  corrections  are  poorly  understood  for  covalent 
semiconductors  like  the  III-V  compounds.  On  the  other  hand,  local  density 
approximation  calculations  g-'<'e  results  v.-hioh  arc  highly  dependent  on  the 
choice  of  the  e>:chc.r.ro-coi'r-;laticn  potential,  .ko  single  potential  has  been 
found  to  give  good  results  for  all  clasrer  of  cempou. is.  For  example, 


Herman,  _e_t  found  that  in  using,  both  the  Slater^  and  the  Kohn-vSnam^* 
exchange  potentials  for  sev:-ra’  serriconduotcra,  the  Slater  potential  agreed 
more  closely  v;ith  experime.-.t  for  II-VI  semiconductors,  while  Kohn-Sharn 
exchange  gave  better  results  for  the  III-V  com.pounds.  Even  in  cases  where 
a  given  potential  yields  the  correct  valence  and  conduction  band 
structures,  sucli  calculations  have  not  al\.'ays  been  able  to  place  tlie  core 
levels  in  their  correct  locatiorus  relative  to  the  valence  bands.  Finally, 
local  density  calculations  have  not  been  successful  in  obtaining  the 
ioninntion  energies  for  electrons  in  most  com.pounds.  V.'c  believe  that  many 
of  these  difficulties  v.rr  ?.  result  of  using  a  Hartrce  potential  which 
includes  a  sel  f-rep/.il s.ion  ter.’’.,  on  the  assumptio.n  that  the  potential  from 
o'lf  electron  in  an  infinite  crystal  will  be  negligible.  As  has  been  shev.n 
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by  Kunz,  et  this  asrun.'^tion  is  not  correct  in  many  cases  and  is  a 
source  of  sionificar.t  error  for  insulators. 

The  new  method,  v;hich  will  be  referred  to  as  the  Hartree-plus  method, 
uses  the  correct,  Hartree  potential,  i-e.  with  the  self-repulsion  removed, 
along  with  a  local  exchange-correlation  potential.  The  Hartree-plus  method 
was  applied  by  Kunz,  ct  al. ^ ,  to  the  solid  rare  gases  and  NaCl  yielding 
bands  which  are  in  good  agreement  with  both  experiment  and  previous  H-F 
plus  correlation  calculations.  The  III-V  semiconductors  provide  a  good 
test  of  the  range  of  applicability  for  the  Hartree-plus  method  since  they 
are  wide  band  semiconductors,  as  opposed  to  the  rare  gases  and  NaCl  which 
arc  narrov;  band  insulators. 

In  Sec. 2,  H-F  theory  is  developed  and  correlation  is  discussed.  The 
Hartree-plus  theory  and  details  of  the  calculations  are  presented  in  Sec. 3 
and  Sec.Jl  re sp-;ctivcly.  In  Sec. 5,  the  results  of  Hartree-plus  band 
calculations  for  AlP,  GaP,  and  Ca/.s  ai'c  presented  and  compared  to 
experiment  and  previous  theoretical  calculations. 

2.  BASIC  THFCh" 


The  systc"!  cf  interest  contains  n  electrons  and  N  nuclei.  Using  the 
Born-Oppenh.cimer  approximation^,  and  neglecting  relativistic  effects,  the 


Hamiltonian  is: 


H'"-'  =  ^  j  Z 

I  4^^-  • 


y'/ 


f/j 


v;herc: 


^  -  -  V/ 


St 


Ttio  energy  is  in  Hydbergs,  c^=2,  uppercase  letters  refer  to  nuclei, 
lowercase  letters  denote  electrons,  and  the  prime  on  the  sun  indicates  that 


the  Eclf-rcpulsion  term,  i=j,  is  excluded.  The  problem  now  is  to  solve  the 


Pape  ^  i 


Schrodinger  equation: 

where  '^'''  is  the  exact  n-elcctron  wavefiir.ction  for  the  excited  state, 

and  the  iTj  are  taken  to  include  both  space  and  spin  coordinates. 

To  simplify  this  problem,  we  approximate  the  exact  v.’avefunction 
with  a  single  Slater  determinant  formed  from  a  set  of  one-electrcn 
orbitals,  (^^(xj)),  i=1,**-,n  and  j=1,***,n: 

The  expectation  value  of  H'"*',  for  the  trial  function  ,  will  be  a 

rigorous  upper  bound  to  the  exact  ground  state  energy  if  the  resulting 


energy  is  stable  against  variation  in  the  jij.  Performing  such  a 
variation,  subject  to  the  constraint  that  the  be  crthonormal ,  yields  an 
equation  for  the 

where  F  is  the  Fcck  operator: 


^  (  A ,  VT,  >  ^  <P^Cx,.) (jf  (/j } 


is) 


In  cq.(3)(  P(2,l)  is  an  exchange  operator,  the  nay  be  used  as 

operators,  and  is  the  Fock-Dirac  density  matrix.  Tlie  second  tern 

is  knov;n  as  the  direct  or  ilartree  term,  and  the  last  term  is  the  exchange 
term.  f.’otc  that  the  prime  is  no  longer  necessary  since  the  Hartree 
self-energy  is  exactly  cancelled  by  the  self-exchange.  Diagonalizing  the 
matrix  Aii  gives  the  standard  liartrec-Fock  equation: 

cF  p-ct:,}  i^) 

The  eigenvalues  of  ec.('>)  are  given  meaning  by  Koopnans'  theorem.^. 
Labelling  occupied  states  v.ith  and  virtu.il  states  with  a=n+1,’’*, 
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the  eicenvalues  arc: 

a  _  ^r^)o _  ^ 

—  t’p 

WO  C’fi 


Csrj> 


where  is  the  energy  of  the  single  Slater  determinant  approximation  to 
the  total  m-elcctron  wavefunction  obtained  by  adding  or  removing  the 
orbital  to  the  n  orbitals  used  to  form  the  H-F  ground  state.  Thus,  and 
6.  are  the  energies  that  an  electron  would  have  if  its  instantaneous 
location  were  independent  of  the  instantaneous  locations  of  the  other 
electrons,  and  if  no  relaxation  were  allowed.  To  improve  on  the  H-F 
eigenvalues  while  retaining  the  independent  orbitals,  one  defines 
correlated  energies  by  replacing  the  approximate  energies  in  eq.(5)  with 
the  exact  energies: 

e;  ^  (o 


w«  *-.‘7  "^tZo 

The  energies  and  are  qu-:.siparti  cle  excitation  energies.  The 
quasiparticles  associated  with  these  energies  ai’e  referred  to  as  holes  and 
conduction  electrons  respectively. 

There  are  tv.’o  basic  approaches  to  finding  the  quasiparticle  energies. 
The  better  understood  of  the  two  is  the  H-F  plus  correlation  method,  which 
is  essentially  a  perturbative  approach.  In  this  method,  eq.(i|)  is  solved 
first,  and  then  the  correlation  effects  are  added  in  as  a  small  number  of 
corrections.  Tnis  approach  to  obtaining  correlated  bands  was  reviewed  and 
developed  in  a  formal  v;ay  by  Pantelidcs,  et  The  correlation 
corrections  are  fairly  well  understood  for  both  metals  and  narrow-band 
wide-gap  materials.  Unfortunately  the  III-V  compounds  do  not  fall  into 
either  of  these  classes.  The  other  approach  is  to  use  an  equation, 
analogous  to  the  H-F  equation,  which  somehow  incorporates  correlation  into 
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the  potential.  Tnis  is  the  approach  taken  by  the  JIartree-plus  nc-thocJ. 

3.  THE  HARTHEE-PLUS  METHOD 

In  the  past,  band  calculations  usin^  local  exchange-correlation 
potentials  have  solved  an  equation  of  the  form: 

Os.r.  (X,)  (pj 


v;here: 

is  a  single  particle  operator  which  replaces  the  nonlocal  exchange  in  the 
Foak  operator  with  a  function,  ,  of  the  local  density ,^y?(r) .  Unlike 

the  Fock  operator,  eq.(7)  no  longer  has  an  exact  cancellation  of  the 
self-repulsion  by  a  self-exchange.  It  has  generally  been  assumed  that  the 
self-repulsion  is  negligible  since  it  comes  from  one  electron  in  an 
infinite  crystal.  This  is  in  fact  true  for  an  electron  in  a  totally 
delocalized  Flcch  state.  However,  v.r.at  is  actually  desired  is  the  energy 
required  to  add  one  quasiparticle  to  the  ground  state  system.  For  all 
insulators,  and  many  semiconductors,  these  quasiparticle  orbitals  are 
local.  In  this  case  it  is  reasonable  to  work  in  the  local  Heitler-Lcndcn 
representation.  For  such  local  orbitals,  the  self-energy  is  not  in  fact 
small  and  should  be  removed  explicitly. 

In  the  local  representation,  let  ^v<(r)  be  the  local  orbital  at 

site  Tnc  corrected  potential  then  is: 


where: 


Vj  Ck)  -  pc  ~~ 


(^) 


is  the  nonlocal  Hartree  potential.  Although  the  Hartrec  potential  is 
nonlocal,  the  nonlocality  is  easy  to  deal  with  by  using  a  single  Hartrce 
potential  for  all  orbitals  of  the  form: 
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where: 

Z.J  Ls.a^  -  7. 

Here  the  bra-ket  notation  is  used  to  emphasize  the  operator  nature  of  the 

s 

self-repulsion  term. 

For  we  have  used  the  local  exchange  potential  of  Kohn-Shan**  and 
the  dielectric  screening  function  of  Robinson,  Dassani,  Knox,  and 
SohriefrerCRB:<S)S.  Thus: 

Y,.  W  is  o) 

where: 

fM  =  I  -  -'(ir)  V/5  Ot--^) 

and: 

cc~.  6146 

Although  it  might  seem  inconsistent  not  to  remove  the  self-exchange  in 
eq.dD,  one  can  argue  that  the  self-exchange  will  only  be  a  significant 
fraction  of  the  total  exchange  potential  in  regions  of  low  density  where 
F(fC)  approaches  zero. 

THE  CALCULATIONS 

To  calculate  the  Hartree-plus  bands  shown  in  Figs.  1-3,  we  first 
performed  an  LCAO  calculation  using  H-F  atomic  wavefunctions  obtained  by 
Bagus,  Gilbert,  and  Roothaan^.  In  this  calculation  Nesbet's  symmetry  and 
equivalence  restrictions^®  were  applied.  Also,  the  crystal  potential  was 
approximated  by  the  first  term  of  its  expansion  in  Yi^t  i-e.  1=0.  Tne 
next  tern  in  the  potential  would  be  1=3  and  should  not  have  a  large  effect 
on  the  s,  p,  and  d  orbitals,  v;hich  are  all  that  are  of  interest  here.  The 
potential  and  matrix  elements  are  evaluated  using  10  shells  about  the 
central  site.  V.'e  represent  the  remainder  of  the  crystal  by  a  residual 
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Kadclung  potential.  To  obtain  this  potential,  vre  used  effective  charges 
given  by  Phillips^^.  This  LCAO  calculation  can  not  be  expected  to  give 
good  results  for  the  conduction  levels  due  to  the  use  of  Kesbet's 
restrictions^®.  The  self-energies  of  the  atonic  orbitals  have  been  removed 
for  both  the  occupied  and  the  unoccupied  parts  of  the  orbitals.  This  means 
that  the  conduction  energies  have  been  obtained  as  though  the  conduction 
electrons  felt  the  potential  from  n-1  electrons,  instead  of  all  n  valence 
electrons.  To  rcnedy  this  the  self-energies  are  explicitly  reinserted  for 
the  virtual  states.  At  this  point  the  Hartree-plus  matrix  is  no  longer 
diagonal.  Hov;evcr  the  coupling  between  the  occupied  and  virtual  states  is 
very  v;eak  and  a  second  diagcnalization  produces  new  valence  bands  which, 
for  GaP,  differ  in  energy  from  the  previous  bands  by  at  most  ^5.  Taking 
all  this  into  consideration,  we  believe  that  the  final  hartree-plus 
operator  so  obtained  is  consistent  with  the  configuration  used.  Finally  we 
form  a  Kartree-plus  matrix  using  atcmic  orbitals  for  the  core  states  and  89 
planewaves  for  the  conduction  and  valence  states,  and  diagonalize  it  to 
obtain  the  final  band  structure.  For  AlP,  v;e  repeated  this  calculation 
using  27,  51,  and  65  planev.-aves. 

5.  RESULTS 

In  all  the  band  calculations  performed  here,  energy  levels  were 
obtained  at  21  nonequivalent  points  in  the  Drillouin  zone  for  an  f.c.o. 
lattice.  The  bands  are  drawn  along  axes  of  high  sy-mnetry  which  connect  the 
points  of  high  sym.’netry  in  the  Brillouin  zone.  The  symmetry  labels  used 
are  consistent  v;ith  Parnenter^^,  In  all  oases  the  origin  has  been  chosen 
to  be  the  cation  site.  In  addition,  the  orbitals  and  their  energies  are 
superscripted  to  indicate  whether  they  are  valence  (v)  or  conduction  (c) 


states. 
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The  valence  band  str-actures  for  AlP^  obtained  using  27,  51,  65,  and  89 
planewaveSy  are  summarized  in  Table  1.  No  experimental  or  theoretical 
values  v;ere  available  for  comparison.  The  bands  obtained  with  £9 
planewaves  are  shoivn  in  Fig.1.  By  comparing  the  results  for  several 
numbers  of  planev;aves,  we  can  estimate  the  size  of  the  errors  due  to 
incomplete  convergence.  For  the  top  two  p-derived  bands,  our  results 
appear  to  be  converged  to  within  0.1  eV.  The  third  p-like  band  is 
converged  to  within  about  0.3  eV,  and  the  separation  of  the  s-band  and  the 
p-bands  is  converged  to  within  roughly  O.^i  cV.  On  the  otherhand,  the  basic 
structure  is  stable,  as  is  the  v/idth  of  the  s-band.  We  assume  that  this 
convergence  information  is  also  valid  for  our  other  calculations. 

In  Tables  2a  and  2b,  the  valence  band  structures  for  GaP  and  GaAs 
shown  in  Figs. 2  and  3  are  compared  to  experiments  and  previous  theoretical 
calculations.  It  is  particularly  interesting  to  compare  the  Hartree-plus 
results  v/ith  the  tv;o  first  prir.ciples  OPW  calculations.  For  both  GaP  and 
GaAs,  the  Hartrcc-plus  method  gives  valence  band  v:idths  which  are  in  better 
agreement  with  experiment  than  those  obtained  in  the  first  principles  OPV.’ 
calculations,  and  are  comparable  in  quality  to  the  empirically  adjusted 
calculations,  especially  when  the  lack  of  convergence  is  taken  into 
consideration.  All  of  the  other  calculations  give  better  results  than  this 
v.'ork  for  the  separation  of  the  s-band  from  the  p-bands.  This  error  in  the 
Hartree-plus  results,  1.0  to  2.0  eV,  is  probably  due  to  incomplete 
convergence  and  relaxation  effects  that  have  not  been  included.  None  of 
the  other  calculations  listed  give  the  position  of  the  d-bands  relative  to 
the  top  of  the  valence  bands.  For  the  Hartree-plus  results,  we  believe 
that  the  errors  in  the  d-bands,  which  are  less  than  15!^,  are  a  result  of 
relaxation  effects,  which  would  be  expected  to  shift  the  d-bands  up 
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relative  to  the  top  of  the  valence  bands. 

In  Table  2c,  the  ionization  enercies  for  the  cation  d-levels  of  GaP 
and  GaAs  are  compared  to  experimental  values.  This  inforr.ation  is  not 
given  by  any  of  the  previous  band  calculations  ncntionec’.  The  differences 
between  theory  and  experir.ent  are  easily  accounted  for  by  relaxation  and 
surface  polarization  effects. 

In  Table  3,  the  lowest  transition  energies  are  compared  to  experiment. 
Clearly  the  Hartree-plus  conduction  bands  are  in  very  poor  agreement  v.'ith 
exporirent.  To  ter.t  for  errors  in  the  treatment  of  the  self-energies,  v;e 
performed  a  calculation  for  AlP  with  the  self-energies  set  equal  to  zero. 
Tnis  calculation  yielded  conduction  bards  v.-hich  differ  from  those  with 
nonzero  sel f-energy  by  no  none  than  0.2  c;V.  ’.v’e  therero>”c;  conclude  that  the 
errors  in  the  conduction  bands  arc  due  to  the  choice  of  the 
excharge-ocrrelation  potential,  It  is  not  really  surprising  that  cur 
cr.cice  of  exor.angc-corrc]  atior-  polcntla  i.  gives  good  renults  for  the  valence 
bands  but  not  foi'  the  conduction  bards.  As  was  shown  by  Pantelides, 

for  both  insulators  and  vnide  gap  .•  tnicenduotors ,  correlation  effects 
in  valence  bands  are  a  result  of  virtual  scattering  of  holes,  while 
correlation  of  the  ccnducticn  bands  comes  primarily  frnrr.  virtual  scattering 
of  electrons.  Since  the  mechanisms  for  correlation  in  valence  and 
conduction  bands  are  different,  it  is  unreasonable  to  expect  one 
excliar.ge-correlation  potential  to  be  valid  for  both  cases.  V.'o  suggest  tiiat 
in  future  ilartrce-plus  calculations,  it  should  be  possible  to  obtain  better 
results  by  using  two  separate  potentials  for  the  valence  and  conduction 
states. 

6 .  CO.'.’ti  I'S.TON 

Pased  upon  the  results  given  by  Kur.z,  et  al.\  for  the  rare  gases  and 


,  for  the  rare  gases 


the  results  for  III-V  semiconductors  reported  here,  v;e  conclude  that  the 
Hartree-plus  method  should  give  good  results  for  the  valence  bands  of  both 
insulators  and  semiconductors  v;hen  used  v/ith  the  RBKS  screened  Kohn-Sham 
exchange.  This  method  not  only  gives  good  results  for  the  valence  band 
widths,  it  also  gives  reasonable  values  for  the  core  d  levels  and  the 
ionization  potentials.  Based  upon  the  calculations  given,  we  suggest  that 
it  might  be  possible  to  find  a  single  exchange-correlation  potential  that 
could  be  used  to  obtain  good  conduction  bands  for  insulators  and 
semiconductors  in  a  two-potential  band  calculation.  Such  a  calculation 
v;ould  use  one  potential  for  valence  states  and  a  second  potential  for 
conduction  states.  This  type  of  calculation  not  only  would  be  expected  to 
give  better  Hartree-plus  band  structures,  but  would  also  be  more  consistent 
with  what  we  know  about  correlation  in  semiconductors  and  insulators. 
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CHAPTER  1 


INTRODUCTION 


The  electronic  structure  of  surfaces  is  important  in  many 
scientific  and  technological  fields,  such  as  corrosion,  semiconductor 
devices,  and  catalysis.  However,  the  understanding  of  fundamental  solid 
surface  properties  has  lagged  behind  the  progress  in  understanding  of 
bul'K  properties.  This  situation  has  been  a  result  of  the  lack  of 
experimental  methods  for  the  quantitative  characterization  of  the 
surface  structure,  and  the  theoretical  complexities  resulting  from  both 
the  loss  of  symmetry  at  the  surface  and  any  possible  reconstruction 
occuring  at  the  surface.  Recently  developed  experimental  techniques 
such  as  electron  energy  loss  spectroscopy  (ELS)  (  ),  two  photon 
picosecond  spectroscopy  (  ),  surface  extended  x-ray  absorption  fine 
structure  (EXAFS)  (  ),  and  extended  appearance  potential  fine 
structure  (EAPFS)  (  )  have  made  accessible  the  experimental  study  of 
the  structural  and  electronic  properties  of  real  surfaces.  The 


( 
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development  of  powerful  computational  methods  has  made  possible  the 
theoretical  study  of  these  systems. 


1.1  THE  INITIAL  ADSORPTION  OF  OXYGEN  ONTO  THE  ALUMINUM  (100)  SURFACE 

Aluminum  surfaces  oxidize  readily  upon  exposure  to  molecular 
oxygen.  Several  experimental  studies  have  been  performed  on  the  Initial 
oxidation  of  aluminum  surfaces  (  ).  Ultraviolet  photoemission 
spectroscopy  (UPS)  (  )  and  Auger  spectroscopy  (  )  studies  have  been 
reported  for  many  crystal  orientations  of  aluminum,  including  the  (100) 
surface.  LEED  (  ),  EXAFS  (  ),  and  EAPFS  (  )  studies  have  been  reported. 

Ultraviolet  photoemission  experiments  (  )  have  reported  a  valence 
band  resonance  at  1.5  eV  below  the  Fermi  level  for  low  oxygen  exposure. 
This  resonance  is  attributed  to  the  oxygen  2p  band.  The  shift  of  the 
aluminum  2p  core  level  has  also  been  measured  for  both  low  and  high 
oxygen  exposures.  This  shift  toward  lower  energy  is  due  to  the  dipole 
moment  arising  from  the  charge  transfer  from  the  aluminum  substrate  to 
the  adsorbed  oxygen.  V/ork  function  measurements  have  also  been  reported 
(  ),  with  the  work  function  found  to  decrease  with  increasing  oxygen 
exposure.  This  decrease  has  been  interpreted  as  being  a  result  of 
incorporation  of  the  electronegative  oxygen  adatoms  beneath  the  aluminum 


surface. 
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LEED,  EXAFS,  and  LAPPS  techniques  are  able  to  probe  the  surface 
geometry  of  the  alur.inux  plus  oxygen  adsorbate  systetn,  and  thus  these 
studies  are  of  special  interest  as  they  will  help  to  differentiate 
between  various  theoretical  models.  For  low  oxygen  coverage,  the 

f* 

reported  Al-0  bond  length  is  1.98  A  (  ),  in  agreement  v;ith  previous 
theoretical  calculations  (  ).  For  higher  exposures,  a  value  of  1.88  A 
is  reported  (  ) ,  which  correspondsto  the  nearest~neighbor  distance  in 
AI2O3 . 

For  theoretical  studies  the  (100)  surface  is  ideal:  calculated 
electronic  properties  can  be  compared  to  existing  Auger  and  UPS  results, 
and  calculated  bond  positions  and  lengths  can  be  compared  with  LEED, 
EXAFS,  and  EAPF3  results. 

Lang  and  Vi'illiams  (4)  have  studied  theoretically  the  adsorption  of 
oxygen  onto  a  jellium  surface.  While  jelliun  is  a  good  approximation  to 
bulk  aluminum,  it  is  not  adequate  to  describe  the  local  adsorption 
process.  For  this  reason,  cluster  studies,  which  are  particularly  well 
suited  to  describe  such  local  phenomena  (  ),  have  been  performed. 
Harris  and  Painter  (  )  and  Kessmer  and  Salahub  (  ) ,  have  studied  the 
adsorption  of  atomic  oxygen  onto  the  hole  site  of  the  A1  (100)  surface, 
reporting  electronic  structures  which  agree  with  existing  UPS  data  for 
the  initial  adsorption  of  oxygen  onto  this  alimunum  surface.  Studies  of 
the  adsorption  of  oxygen  at  the  remaining  two  high-symmetry  sites  of 
this  surface  have  not  previously  been  reported.  An  accepted  explanation 
of  the  incorporation  of  oxygen  atoms  beneath  the  aluminum  surface  does 
not  yet  exist:  the  calculated  electronic  structure  and  binding 


potential  curves  for  the  on-top  and  bridge  sites  presented  here  should 
shed  light  on  this  problem. 

In  this  study,  the  interaction  between  adsorbate  atoms  (in  this 
case  oxygen)  and  the  aluminum  (100)  surface  is  theoretically 
investigated  by  modeling  the  surface  with  a  small  number  of  atoms  and 
then  using  the  unrestricted  Hartree-Fock  approximation  to  calculate  the 
electronic  structure  and  potential  energy  curves.  From  these  potential 
energy  curves  the  binding  energy,  bond  distance,  and  force  constants  of 
the  interaction  are  determined.  The  calculated  one-electron  orbitals 
give  added  insight  into  the  chemisorption  process,  and  calculated 
densities  of  states  are  compared  with  photoeraission  data.  This 
calculational  procedure  is  able  to  accurately  describe  the  localized 
bonding  of  an  adsorbate  onto  a  surface  site,  the  aspect  of  chemisorption 
considered  most  important  (  ). 

The  localized  nature  of  the  chemisorptive  bond  is  indicated  by 
several  kinds  of  experimental  evidence.  An  adsorbed  atom  or  molecule 
can  be  observed  to  hop  from  one  localized  site  to  another  (  ). 
Infrared  spectroscopy  studies  have  found  that  the  vibrational  spectra  of 
Intermediates  on  surfaces  are  often  very  similar  to  the  spectra  of 
isolated  molecules  (  ) .  These  experimental  results  reinforce  the 
intuitive  notion  that  the  chemisorptive  bond  is  similar  in  nature  to  the 
familiar,  localized  chemical  bond. 
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In  general,  the  loss  of  symmetry  in  the  direction  perpendicular  to 
the  surface  leads  to  the  choice  between  two  drastic  approximations  in 
theoretically  modeling  the  surface-adsorbate  interaction.  The 
maintenance  of  translational  symmetry  in  the  plane  parallel  to  the 
surface  greatly  simplifies  the  computational  difficulties  and  allows  for 
the  use  of  the  powerful  methods  developed  for  bulk  band  structure 
calculations.  The  alternative  approximation,  simulating  a  small  section 
of  the  surface  by  a  cluster  of  atoms,  can  more  effectively  describe  the 
localized  interaction  of  the  surface  with  a  single  atom. 

Numerous  calculations  of  the  electronic  structure  of  surfaces  using 
the  approximation  of  a  semi-infinite  surface  have  been  reported  (  ). 
The  major  drawback  of  this  approach  is  that,  in  general,  a  semi-infinite 
monolayer  of  adsorbed  atoms  must  be  considered  in  a  chemisorption 
calculation.  This  lim.itation  is  a  direct  result  of  the  symmetry 
assumed,  and  therefore  only  qualitative  calculations  can  treat  a  single 
adatom  on  a  semi-infinite  surface. 

By  relaxing  the  symjnetry  requirements  of  the  semi-infinite  surface, 
one  gains  increased  flexibility  and  freedom.  For  instance,  atoms  can  be 
brought  down  over  a  variety  of  possible  bonding  sites.  Roughened  and 
stepped  surfaces,  exposed  corners,  and  small  particles,  which 
experimental  work  indicates  are  often  sites  favored  for  chemisorption, 
can  easily  be  simulated  by  the  cluster  approach.  Use  of  the 
semi-infinite  surface,  by  contrast,  is  limited  to  studies  of  the  perfect 


surface 
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Once  the  appropriate  cluster  has  been  chosen  to  simulate  the 
surface,  the  electronic  structure  of  the  cluster-adsorbate  system  must 
be  calculated.  The  method  used  here  is  the  unrestricted  Hartree-Fock 
approximation,  from  v/hich  electronic  wavef unctions,  electron  energy 
eigenfunctions,  and  the  total  electronic  energy  are  calculated.  This 
calculation  is  repeated  at  varying  adsorbate-surface  separations,  and 
the  potential  energy  curve  can  be  generated  from  the  differences  in 
total  electronic  energy. 

The  Hartree-Fock  approximation  as  the  method  of  calculation  is 
attractive  due  to  two  important  advantages.  The  method  is  a  convenient 
first  step  towards  an  £b  initio  solution  of  the  mv^ny-particle 
Schrodinger  equation  for  a  fernion  system,  and  gives  a  useful 
zeroth-order  vtavefunction  for  a  perturbation  calculation  of  the 
many-body  problem.  Second,  it  is  the  most  advanced  calculational  method 
in  which  elementary  physical  intuition  is  applicable.  In  more  complete 
calculations ,  there  is  no  longer  a  one-to-one  correspondence  between 
particles  and  one-particle  wavefunctions,  and  the  independent  particle 
approximation  is  no  longer  applicable. 


1.2  THE  ELECTKONIC  PROPEHTIES  OF  THE  SILICA  SURFACE 


The  electronic  structure  of  crystalline  silica  (Si02)  is  of 
technological  interest  in,  for  example,  the  manufacture  of  solid  state 
electronic  devices.  Consequently,  there  has  been  considerable 
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experimental  effort  in  this  area.  In  a  recent  review  of  experimental 
and  theoretical  results,  Griscom  (  )  has  compiled  much  of  the  existing 
data  on  the  bulk  properties  of  silica.  The  band  gap  has  been  found  to 
be  about  8.9  eV,  with  a  valence  band  width  of  approximately  11  eV. 
Ultra-violet  photo  emission  spectroscopy  (UPS)  studies  by  Ibach  and  Rowe 
(  )  have  detected  no  occupied  surface  states  v/ithin  the  band  gap. 
Williams  (  )  has  studied  silica  surfaces  grown  on  crystalline  Si  and 
found  traps  at  2  eV  below  the  bulk  silica  conduction  band  edge.  These 
traps  act  like  Coulomb  centers  with  a  positive  charge.  It  is  with  the 
silica  surface  and  its  defects  that  the  current  investigation  is 
concerned. 

Several  recent  experimental  studies  have  investigated  surface 
states  of  silica  (  ).  Schwidtal  (  )  has  found  a  radiation-induced 
feature  at  SI  eV  on  the  high-energy  side  of  the  Si  L23  VV  (V=valence) 
transition  in  the  Auger  electron  spectrum  (AES),  and  suggests  that  it  is 
due  to  an  Si  L22  VD  (D=defect)  transition.  He  further  suggests  that 
this  defect  is  the  center,  which  is  a  dangling  singly-occupied  sp^ 
orbital  of  silicon.  Ibach  and  Rowe  (  ),  Fujiwara  and  Ogata  (  ), 
Lieske  and  Uezel  (  ) ,  and  Bermudez  and  Ritz  (  )  have  measured  electron 
energy  loss  spectra  (ELS)  and  found  peaks  in  the  second  derivative 
spectrum  at  about  3.5,  5.0,  and  7.4  eV  (Fig.  1).  Ibach  and  Rowe  have 
attributed  these  transitions  to  a  partially  oxidized  surface  region, 
SiOjj,  where  UX'2.  Fujiv;ara  and  Ogata  (  )  have  shown  the  states 
associated  with  these  transitions  to  be  at  the  surface,  and  concluded 
they  were  due  to  metastable  SiOj^.  Lieske  and  Hezol  (  )  have  associated 
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the  final  states  with  Si-Si  bonds  found  in  SiO^^,  Most  recently, 
Bermudez  and  Ritz  (  )  have  studied  this  transition,  and  found  the  final 
state  to  be  due  to  “a  chemically  stable  species  formed  after  rupture  of 
the  silica  network,"  and  attributed  these  properties  to  a  surface  Si=0 
double  bond.  For  this  bond  to  exist,  the  silica  surface  must  undergo 
reconstruction.  Evidence  for  such  reconstruction  from  low  energy 
electron  diffraction  (LEED)  measurements  is  due  to  Janossy  and  Menyhard 
(  ).  Also,  Hochstrasser  and  Antonini  (  )  have  observed  what  they 
believe  to  be  recombination  luminescence  due  to  the  rearranging  of  the 
surface  Si-0  bonds  immediately  after  cleaving.  They  measured  the 
lifetime  of  this  luminescence  to  be  loss  than  lO"^  seconds. 

Theoretical  studies  of  the  silica  bulk  electronic  properties  have 
bean  perfcrned  (  ),  and  the  calculated  densities  of  states  are  in 
good  agreement  with  experiment,  except  for  the  AES  peak  at  91  eV  and  the 
ELS  peaks  at  3.3,  5.0,  and  6.8  eV.  Both  of  these  discrepancies  can  be 
accounted  Vo"  by  allowing  for  the  existence  of  islands  of  silicon  in  the 
bulk  Si02.  Bennett  and  Roth  (  )  have  calculated,  using  the  Huckel 
approximation,  the  electronic  properties  of  many  clusters  which 
approximate  silica,  and  in  particular  gave  attention  to  defects.  These 
calculations  were  performed  in  the  bulk;  however,  the  clusters  were 
sufficiently  small  that  they  can  be  interpreted  in  terms  of  surface 
defects,  or  in  terms  of  a  surface  SiOj^  region.  For  oxygen  deficient 
clusters  v/ith  a  relaxation  of  the  silicon  atoms  near  the  defect,  the 
calculated  energy  levels  showed  transitions  at  approximately  the  same 
energy  values.  These  are,  however,  probably  due  to  small  clusters  (2-3 


atoms)  of  silicon  in  the  SiOj^.  Laughlin,  et  al.  (  )i  using  a  tight 
binding  method,  have  calculated  similar  defects  and  have  also  found 
states  which  appear  in  the  gap,  at  about  5  eV  above  the  valence  band 
edge.  They  have  also  studied  the  E*  center,  and  found  a  level  lying 
just  below  the  bulk  conduction  band.  Pollmann  and  Pantelides  (  )  have 
performed  a  calculation  using  a  Green's  function  formulation  on  the 
ideal  terminated  cubic  ^-cristobalite)  surface.  Using  an  admittedly 
crude  model  for  the  surface,  they  found  no  states  in  the  band  gap. 

The  purpose  of  the  current  investigation  is  to  study  the  silica 
surface,  giving  particular  attention  to  the  three  ELS  peaks  observed  at 
3.5,  5,  and  7  eV.  The  cluster  model  is  again  employed,  for  the  same 
reasons  given  in  section  1.1.  Of  the  three  models  proposed  by 
experimentalists  to  account  for  these  transitions,  the  only  one  which 
has  previously  been  calculated  theoretically  is  theSiOj^  model.  In 
chapter  ^  are  presented  calculations  of  the  one-electron  energy  levels 
for  the  remaining  two  models,  the  E’  center,  and  the  Si=0  surface 
double  bond.  A  comparison  is  then  made  between  the  three  models  and 


experiment. 


METHOD  OF  CALCULATION 


An  exact  solution  to  the  wave  equation  for  a  large  cluster  of  atoms 
is  not  at  present  attainable.  In  fact,  such  a  solution,  while  nice  in  a 
mathematical  sense,  would  be  difficult  to  interpret  physically.  The  ab 
initio  Hartree-Pock  theory  provides  an  approximate  method  of  solution  to 
determine  the  electronic  structure  of  a  large  cluster  of  atoms  which  is 
both  easily  handled  by  contemporary  computers  and  well  suited  for 
straightforward  interpretation. 

2.1  HARTREE-FOCK  THEORY 

Ab  IniijUi  Hartree-Fock  theory  is  used  to  calculate  the  electronic 
wavefunotions  and  energies  for  the  finite  clusters  studied  in  this  work. 
The  solids  studied  contain  atoms  of  low  atomic  number;  therefore 
relativistic  effects  are  small  and  will  be  ignored.  Using  the 
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nonrelativistic  Schrodinger  theory,  the  many-body  wavefunction  T,  which 
is  a  function  of  electronic  and  nuclear  space-spin  coordinates,  is 
determined  by  the  equation 


where  H  is  the  many-body  Hamiltonian 


J 
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Here  we  have  used  atomic  units,  where^h"^,  mg,  and  e=1,  the  unit  of  length 
is  the  bohr  (1  bohr;:^  0.529A) ,  and  the  energy  is  given  in  Hartree  (1  Hy 
'iC27.2  eV).  Upper  case  letters  refer  to  nuclear  properties:  "Rj  is  the 
position  of  the  nucleus,  and  Mj  is  its  mass.  Lower  case  letters 
refer  to  electronic  properties;  'r^  is  the  coordinate  of  the  i^*^ 
electron.  The  four- vector  denotes  both  the  spatial  coordinates  and 
the  spin  of  the  i^h  electron. 

The  Schrodinger  equation  is  simplified  by  employing  the 
Born-Oppenheimer  approximation  (  )  to  separate  the  nuclear  and 
electronic  coordinates.  The  total  wavefunction  is  assumed  to  be 
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separable: 

i  o,h'-  (aA) 

The  nuclear  wavefunction  '{'(X),  is  now  a  solution  of  the  equation 


while  the  electronic  wavefunction, '^(x,X) ,  is  a  solution  of  the  equation 

A  solution  of  this  electronic  Schrodinger  equation  depends  on  the 
nuclear  coordinates  R  and  in  turn  defines  a  potential  energy  function  of 
the  coordinates: 


which  determines  the  motion  of  the  nuclei.  The  problem  is  now  reduced 
to  solving  the  electronic  Schrodinger  equation  for  a  particular  set  of 
nuclear  coordinates.  This  equation  cannot  in  general  be  solved  exactly; 
one  therefore  takes  refuge  in  the  independent  particle  model  for  a 
method  of  obtaining  an  approximate  solution. 

This  approach  is  mathematically  equivalent  to  a  separation  of 
variables  technique.  The  many-electron  wavefunction,  which  is  a 
functional  of  the  space-spin  coordinates  of  all  of  the  electrons,  is 
assumed  to  be  a  product  of  orbitals  which  are  functions  of  one-electron 
space-spin  coordinates.  However,  the  Pauli  principle  requires  that  the 
many-electron  wavefunction  be  antisymmetric  under  the  exchange  of  any 
two  electrons.  Therefore,  the  total  electronic  wavefunction  is 
approximated  by  an  antisyrametrized  product  of  one-electron  orbitals: 
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where  ^is  the  antisynunetrizer.  Since  the  Schrodinger  equation  is  still 
not  in  general  solvable  analytically,  the  variational  principle  is 
introduced  to  obtain  approximate  solutions.  The  variational  principle 
states  that  the  normalized  expectation  value  of  the  Hamiltonian  using 
approximate  wavefunctions  is  a  rigorous  upper  bound  to  the  exact  energy. 
The  approximate  wavefunction  can  then  be  varied  to  minimize  the  energy 
and  thereby  be  the  best  approximation  to  the  exact  wavefunction. 
Applying  the  variational  principle  to  the  electronic  Schrodinger 
equation,  and  using  an  antisymmetrized  product  of  one-electron  orbitals, 
we  obtain  the  Hartree-Fook  equations: 


where 


i;. 

1 


the  Fock  operator: 


The  matrix  elements  are  Lagrange  multipliers,  and  have  been 
introduced  to  insure  the  orthonormality  of  the  one-electron  orbitals. 
In  general,  a  unitary  transformation  can  be  performed  on  the  orbital 
space  which  diagonalizes  this  matrix.  The  Hartree-Fock  equations  can 
then  be  written  in  the  standard  form: 


These  integro-differential  equations  are  normally  salved  via  an 
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iterative  procedure.  Thus,  the  electron  orbitals  are  determined  in  the 
average  field  of  all  the  other  electrons.  The  only  constraints  imposed 
on  the  wavefunction  are  that  the  orbitals  be  orthonormal  and  be 
functions  of  the  space-spin  coordinates  of  only  one  electron.  This 
level  of  approximation  is  called  generalized  Hartree-Fock  (GHF)  and  has 
yet  to  be  solved,  so  further  constraints  are  imposed  to  simplify  the 
calculation.  The  simplest  constraint  which  can  be  imposed  is  to  require 
that  the  orbitals  be  eigenfunctions  of  S^,  i.e.,  the  orbitals  be  of  the 
form: 


't  ^  )  f 


where  ‘‘(r)  is  the  spatial  part  of  the  orbital,  and  u  and  ^  are  the 
spin-up  and  spin-down  eigenfunctions  of  S^,  respectively.  The  form  of 
the  Hartree-Fock  equations  is  unaffected  by  this  constraint,  and  this 
level  of  approximation  is  known  as  the  unrestricted  Hartree-Fock  (UHF). 
The  UHF  method  has  been  used  for  numerous  practical  calculations  of 
atoms,  molecules,  and  solids  (  ),  and  is  the  method  used  in  this 
study . 

Additional  restrictions  may  be  placed  upon  the  one-electron 
orbitals  in  order  to  further  simplify  the  calculations.  The  spatial 
parts  of  the  orbitals  may  be  required  to  be  symmetry  eigenfunctions  of 
the  cluster  under  consideration,  i.e.,  each  orbital  must  transform  as  an 
irreducible  representation  of  the  space  group  of  the  cluster.  Finally, 
for  doubly  occupied  orbitals  the  spin-up  and  spin-down  spatial  orbitals 
may  be  assumed  to  be  equivalent.  This  level  of  approximation  is  called 
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the  restricted  Hartree-Fock  method  (RHF).  One  advantage  of  RHF  is  that 
the  one-electron  orbitals  are  eigenfunctions  of  both  spatial  symmetry 
and  spin.  However,  for  singly  occupied  orbitals  complications  arise  in 
that  the  Fock  operator  differs  for  open  and  closed  shells  of  the  same 
symmetry,  thus  not  allowing  the  single  determinant  to  properly  describe 
the  wavefunction.  The  UHF  method,  while  not  necessarily  obtaining  a 
wavefunction  which  is  an  eigenfunction  of  spin  and  spatial  symmetry, 
does  give  a  lower  (i.e.,  better)  eigenfunction  of  energy.  Since  the 
"best"  solution  to  equation  2.6  is  the  wavefunction  which  gives  the 
lowest  energy,  the  UHF  method  has  therefore  been  used  in  this 
calculation. 


2.2  CLUSTER  MODEL 


Localized  phenomena  in  solids,  such  as  point  defects  or  surface 
adsorbates,  can  be  modeled  effectively  by  using  a  finite  cluster  of 
atoms  (  ) .  Since  the  cluster  is  intended  to  represent  a  much  larger 
system,  one  must  apply  appropriate  boundary  conditions  for  the  cluster 
to  simulate  the  effects  of  the  environment.  These  have  been  discussed 
by  Kunz  and  Klein  (  ),  and  will  be  briefly  reviewed  here. 

The  solution  to  the  Hartree-Fock  equation  (2.11)  is  sought  for  a 
finite  cluster.  Suppose  the  region  of  the  cluster  is.  termed  A,  and  the 
remainder  of  the  system,  the  environment  of  A,  is  termed  E.  The  problem 
is  to  partition  the  system  rigorously  into  a  cluster  and  an  environment. 
Let  us  consider  the  method  of  local  orbitals  of  Adams,  Gilbert,  and 
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Kunz(  ).  Here  the  Fock  operator,  F,  of  the  entire  system  is  divided 
into  a  part  which  includes  the  kinetic  energy,  nuclear  attraction  of  the 
electrons  and  nuclei  inside  A,  along  with  the  electron-electron 
potential  for  the  electrons  assigned  to  A,  F^^;  and  another  part  which 
is  the  potential  in  A  due  to  the  environment,  U^: 


It  is  desired  to  study  only  part  of  the  system,  so  instead  of  the  normal 
Hartree-Fock  equation  (2.11),  consider  instead  the  local  orbitals 
equation: 


where  W  is  an  arbitrary  Heraitian  operator.  Consider  now  an  ionic 
system.  The  potential  due  to  the  environment,  U^,  may  be  divided  into 
two  parts:  V‘^‘  is  an  ionic  (Madelung)  contribution  and  is  long  range, 

C  O 

and  is  the  remainder  and  is  short  range.  Let  W  =  V^.  Making  use 
of  the  projector  properties  of  ,  and  solving  for  the  occupied 
orbitals,  equation  (2.1l|)  becomes 


Considering  only  the  orbitals  of  equation  (2.15)  which  lie  in  A,  the 
solutions  should  penetrate  only  weakly  into  E.  Because  does  not 
appreciably  penetrate  A,  and  because  v|  is  cancelled  by  i  in  the 
limit  of  self-consistency,  the  approximate  equation  for  the  cluster, 
including  the  interaction  with  the  remainder  of  the  system,  is  then 
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This  approximate  equation  is  simply  the  equation  for  the  wavefunction  in 
the  cluster  of  interest  imbedded  in  the  Madelung  field  of  the  remainder 
of  the  solid.  The  eigenvalues of  equation  (2.14)  or  the  approximate 
equation  (2.16)  represent  the  eigenvalues  Ci  of  the  infinite  solid  only 
in  the  limit  that  the  orbitals  '|>j^  are  localized  in  the  region  A. 


A  basis  function  expansion  is  used  to  determine  the  unrestricted 
Hartree-Fock  orbitals.  Here  gaussian  type  orbitals  have  been  used  for 
this  expansion.  This  type  of  function  offers  the  advantage  of  having 
analytic  solutions  for  the  required  integrals;  however,  many  such 
functions  are  required  to  achieve  the  desired  accuracy.  Cartesian 
gaussian  functions  have  been  used.  The  spatial  part  of  these  functions 
has  the  form: 


J  H 


J 

The  determination  of  the  gaussian  exponents,  for  atoms  and  ions  has 
been  the  subject  of  considerable  study,  and  tabulated  sets  exist  in  the 
literature.  An  all-electron  calculation  including  the  core  electrons  of 
aluminum  and  silicon  would  be  impractical,  so  the  ab  initio  effective 
potential  of  Topiol,  et  al.  (  )  is  used  to  replace  these  core 
electrons.  The  integrals  are  performed  using  the  standard  POLYATOM 
integrals  program,  and  the  unrestricted  Hartree-Fock  calculations  are 
done  with  the  G.  T.  Surratt  program  UHFONE. 


Page 


THE  INITIAL  ADSORPTION  OF  OXYGEN  ONTO  THE 
ALUMINUM  (100)  SURFACE 


AluniriUn  crystalizes  in  the  FCC  form,  with  a  lattice  constant  of 
11.05  A.  Aluninum  is  a  metal,  and  the  bulk  electronic  structure  is 
approximated  very  well  by  the  "jellium"  model.  However,  attempts  to 
employ  this  model  to  describe  the  local  chemisorption  process  (10)  have 
produced  results  which  are  not  consistent  with  experiment  (6).  Cluster 
techniques  have  been  used  to  study  the  chemisorption  process  (11,12,40), 
with  results  v/hich  agree  well  with  experiment  if  care  is  taken  in 
choosing  the  cluster. 

Chemisorption  is  a  localized  process,  and  as  such  the  cluster  model 
should  provide  an  accurate  description.  The  difficulties  that  one 
encounters  here  are  resultant  from  the  fact  that  aluminum  is  a  metal, 
and  as  such  is  not  well  described  by  a  small  cluster  of  atoms  (fig. 
3.1).  One  must,  therefore,  take  into  account  the  size  of  the  cluster 
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used  to  describe  the  substrate.  These  effects  will  be  discussed,  along 
with  the  binding  energy  of  the  oxygen  atom  for  each  surface  site 
considered.  Calculated  one-electron  energy  levels  for  these  binding 
sites  are  compared  with  the  results  of  photoemission  experimemts. 


3.1  RESULTS  OF  CALCULATIONS 


The  positions  of  the  aluminum  centers  in  the  clusters  are  chosen  to 
reproduce  the  unreconstructed  structure  of  the  aluminum  (100)  surface. 
It  has  been  assumed  here  that  the  lattice  does  not  undergo  signifigant 
distortions  at  the  surface.  Tnis  has  been  shown  to  be  true  from 
comparison  of  experimental  and  theoretical  LEED  studies  of  clean 
aluminum  surfaces  (4^0. 

As  possible  surface  sites  for  the  adsorbed  oxygen  atom,  considered 
here  are  the  three  high-symmetry  points  of  the  (100)  surface  (fig. 
3.2).  These  are  a  position  directly  above  an  aluminum  center  of  the 
first  substrate  layer  (on-top  position);  the  position  central  to  four 
aluminum  atoms  of  the  first  layer,  above  a  second  layer  atom  (hole 
position);  and  above  (or  below)  the  midpoint  of  a  line  connecting  two 
nearest  neighbor  aluminum  centers  (bridge  position). 

The  interaction  of  oxygen  with  the  aluminum  surface  at  the  on-top 
position  has  been  modeled  here  using  two  clusters:  AlO,  with  a  single 
aluminum  atom  representing  the  surface;  and  AlgO,  with  five  first  layer 
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atoms  and  four  second  layer  atoms  modeling  the  surface.  The  binding 
energy  of  the  oxygen  atom  to  the  aluminum  surface  is  shown  as  a  function 
of  its  distance  above  the  surface  in  figure  3.3*  For  both  surface 
clusters,  the  equilibrium  distance  is  about  3-5  bohr.  The  binding 
energies  differ  between  the  two  clusters;  this  energy  is  2.25  eV  for 
the  AlO  cluster,  and  1.36  eV  for  the  AlgO  cluster.  The  charge  transfer 
from  the  aluminum  surface  to  the  oxygen  atom  is  about  0.7  electron  in 
both  cases. 

For  the  hole  site,  the  aluminum  surface  is  modeled  with  tv/o 
clusters:  AI5O,  which  has  four  aluminum  centers  in  the  first  surface 
layer,  and  the  central  aluminum  in  the  second  layer;  and  AlgO,  which  is 
the  AI5O  cluster,  to  which  four  additional  second  layer  aluminum  centers 
are  added.  The  importance  of  including  the  second  layer  atoms  at  this 
site  has  been  shown  by  previous  theoretical  studies  of  the  interaction 
of  oxygen  with  metal  surfaces  (^(5).  The  oxygen  binding  energies  as  a 
function  of  the  distance  above  the  surface  plane  are  shown  for  these 
clusters  in  figure  3*^»  Both  clusters  have  a  minimum  total  energy 
(maximum  binding  energy)  when  the  oxygen  center  is  about  0.2  bohr  below 
the  plane  of  the  first  surface  layer.  The  cluster  with  only  one  second 
layer  aluminum  atom  is  found  to  bind  the  oxygen  adsorbate  by  5.^  eV, 
with  a  charge  transfer  to  the  oxygen  of  1.^  electrons.  The  larger 
cluster  is  found  to  bind  the  oxygen  by  ^.8  eV,  with  a  charge  transfer  of 
1.3  electrons  to  the  adsorbed  oxygen. 
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The  binding  energies  computed  for  this  site,  as  well  as  those  for 
the  bridge  site,  must  be  carefully  extracted  from  the  calculated  total 
energies.  Since  the  oxygen  atom  penetrates  into  the  cluster  of  aluminum 
centers,  the  basis  functions  used  to  describe  the  electronic 
wavefunctions  associated  with  the  oxygen  adorbate  may  improve  the 
description  of  the  substrate  electronic  structure,  thus  lowering  its 
total  energy.  This  improvement  in  the  substrate  wavefunction  can  be 
integrated  into  the  caculation  of  the  binding  energy  by  simply  including 
the  basis  functions  of  the  oxygen  atom  when  calculating  the  total  energy 
of  the  substrate  cluster. 

The  bridge  site  has  been  modeled  here  with  three  different 
clusters.  The  first  cluster  consists  of  two  nearest  neighbor  aluminum 
centers  of  the  surface,  anc  an  oxygen  atom  above  the  midpoint  of  the 
line  connecting  these  centers.  The  second  cluster  contains  these  two 
centers,  as  well  as  two  additional  centers  beneath  them  from  the  second 
aluminum  layer.  For  the  final  cluster,  tv;o  aluminum  centers  from  the 
third  layer  are  added.  The  inclusion  of  the  deeper  layer  aluminum 
centers  at  this  bonding  site  is  necessary  because  the  oxygen  atom 
penetrates  the  surface,  and  can  move  vertically  in  a  "tunnel*'  between 
pairs  of  aluminum  centers  (see  fig.  3.2).  The  binding  energies  of  an 
oxygen  atom  interacting  vdth  these  aluminum  clusters  are  showTi  in  figure 
3.5.  For  the  cluster  with  two  aluminum  centers,  the  oxygen  adsorbate  is 
bound  1.3  bohr  above  the  plane  of  the  surface.  A  charge  of  1.25 
electrons  is  transfered  to  the  oxygen  atom,  which  is  bound  by  3.1  eV. 
With  the  inclusion  of  the  two  second  layer  aluminum  centers  in  the 
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cluster,  the  equilibrium  position  of  the  oxygen  center  moves  to  1.65 
bohr  beneath  the  surface,  which  is  just  above  halfv/ay  between  the  two 
aluminum  layers.  Tne  binding  energy  is  found  to  be  3.8  eV,  with  a 
charge  transfer  of  1.^  electrons  to  the  oxygen.  Inclusion  of  the  third 
aluminum  layer  produces  no  signifigant  change  in  these  results; 
however,  this  cluster  does  show  that  the  oxygen  center  docs  not  want  to 
move  deeper  into  the  bulk  material,  as  the  total  energy  of  the  system  is 
found  to  be  lowest  for  an  oxygen  center  position  just  below  the  first 
surface  layer. 


3.2  DISCUSSION  OF  RESULTS 


For  the  initial  stages  of  adsorption,  the  oxygen  adsorbate  is 
clearly  bound  below  the  aluminum  surface.  The  binding  is  strongest  at 
the  hole  site.  An  electric  charge  of  about  1.35  electrons  is  also 
transfered  to  the  adso.^bed  oxygen.  These  results  are  in  agreement  with 
previously  reported  work  function  studies,  in  which  a  reported  decrease 
in  the  work  function  of  the  aluminum  surface  with  increasing  oxygen 
exposure  (up  to  one  monolayer  of  coverage)  has  been  attributed  to  the 
penetration  of  negative  oxygen  ions  beneath  the  aluminum  surface  (6). 

Tne  ultra-violet  photoeraission  spectrum  has  been  reported  for  low 
oxygen  coverage  on  the  aluminum  (100)  surface  (5);  this  spectrum  is 
shown  in  figure  3.6.  Also  shovm  is  the  projected  density  of  states  for 
the  oxygen  adsorbate  previously  calculated  by  Messmer  and  Salahub  (12) 


for  an  oxygen  position  at  the  aluminum  surface  at  the  hole  site.  These 
authors  did  not  calculate  the  electronic  structure  for  an  oxygen 
position  below  the  surface,  nor  did  they  study  the  interaction  of  an 
oxygen  adsorbate  with  any  of  the  other  high-symmetry  (100)  surface 
sites. 

As  reported  elsewhere  (MO) ,  the  one-electron  energy  levels  of  this 
calculation  are  in  agreement  with  experiment  and  previous  theory.  These 
energy  levels  are  shown  in  figure  3.6  for  both  the  hole  and  bridge 
sites,  with  the  oxygen  position  corresponding  to  the  minimum  total 
energy  for  the  cluster.  For  both  sites,  the  energy  levels  reported  here 
are  those  of  the  largest  cluster  calculated. 

Photoemission  experiments  studying  the  aluminum  core  levels  at  low 
oxygen  exposures  have  been  reported,  with  a  reported  shift  in  the 
aiuminura  2p  energy  of  1.3eV  tov;ard  lower  energy  (5).  In  the  present 
calculation,  the  shift  of  the  alurainm  core  levels  is  found  to  be  1.1  eV 
in  the  same  direction.  Tnc  cause  of  this  core  level  shift  is  the 
electric  dipole  produced  by  the  charge  transfer  from  the  aluminum 
substrate  to  the  oxygen  atom  upon  adsorption.  It  has  been  noted  by 
Flodstrbm,  e^  al.  that  the  aluminum  core  level  shift  increases  to  a 
value  of  2.6  eV  at  higher  oxygen  exposure  (about  one  monolayer  of 
coverage),  thus  lowering  these  orbital  energies  to  the  values  found  in 
AI2O3.  These  results  are  evidence  for  a  two  stage  oxidation  process, 
with  the  first  step  being  chemisorption  at  the  hole  site,  and  the  second 
step  being  incorporation  of  oxygen  beneath  the  aluminum  surface. 
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Hie  agreement  between  the  calculated  one-electron  energies  for  the 
hole  site  presented  here,  along  with  those  of  Messmer  and  Salahub,  and 
reported  experimental  UPS  spectra  tend  to  support  the  idea  that  the 
initial  chemisorption  of  oxygen  occurs  at  the  hole  site,  as  one  would 
expect  from  the  binding  energies  reported  in  section  3.1. 

The  hole  site  is  also  favored  by  extended  appearance  potential  fine 
structure  (EAPFS)  measurements.  Using  this  technique,  den  Boer,  ^  a^. 
(i4)  have  reported  a  nearest  neighbor  oxygen-aluminum  distance  of  1.98  K, 
or  3.7^  bohr.  Messmer  and  Salahub  have  reported  a  value  of  2.02  K  (3.82 
bohr),  and  tlie  results  of  the  calculations  presented  here  place  the 
oxygen  at  1.92  Jl  (3.63  bohr)  from  the  nearest  aluminum  center  for  the 
hole  site.  For  higher  oxygen  exposure,  den  Boer^et  al^.  have  reported 
an  experimental  value  of  this  nearest  neighbor  distance  of  1.91  K.  This 
figure  is  probably  an  average  of  the  two  nearest  neighbor  spacings  found 
in  AI2O2 ;  these  being  1.85  A  (3.bl  bohr)  and  1.97  ^  (3*72  bohr).  Tne 
longer  of  these  spacings  corresponds  to  oxygen  binding  at  the  hole  site, 
while  the  shorter  corresponds  to  oxygen  adsorption  at  the  slightly  less 
energetically  favorable  bridge  site.  The  present  calculation  places  the 
nearest  neighbor  distance  at  the  bridge  site  at  1.71  A,  or  3.25  bohr. 

Tne  results  presented  above,  along  with  previous  theoretical  and 
experimental  results,  show  that  for  the  initial  interaction  of  oxygen 
with  the  aluminum  (100)  surface,  chemisorption  takes  place  at  the  hole 
site.  Ihe  present  studies  have  extended  previous  theoretical 
understanding  of  the  interaction  of  oxygen  with  this  surface  to  include 
the  two  additional  high-synetry  (100)  sites.  Uork  function,  UPS,  and 
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EAPFS  studies  have  reported  that  the  nature  of  this  interaction  changes 
when  the  oxygen  coverage  approaches  one  nonolayer;  the  theoretical 
results  presented  here  lead  one  to  conclude  that  at  this  coverage, 
oxygen  begins  to  adsorb  at  the  bridge  site,  since  the  energitically  more 
favorable  hole  sites  are  already  occupied.  The  chemisorption  of  oxygen 
at  the  bridge  site  allows  for  the  incorporation  of  the  oxygen  adatoms 
beneath  the  aluminum  surface  at  high  levels  of  oxygen  exposure,  leading 
to  the  formation  of  the  oxide  AI2O3. 


THE  ELECRONIC  PROPERTIES  OF  THE  SILICA  SURFACE 


Silica  occurs  in  six.  different  cry.stalline  forms,  as  well  as  the 
amorphous  form.  All  forms  of  Si02  are  based  on  the  SiOtj  tetrahedron, 
and  differ  in  the  Si-O-Si  bonding  angle.  The  crystalline  form  most 
often  studied  experimentally  is  pC-quartz,  and  it  is  this  form  which  we 
have  used  here.  <j(-quartz  is  hexagonal  in  structure,  with  three  Si02 
molecules  in  each  unit  cell.  The  Si-O-Si  bond  angle  is  . 

The  silica  surface  may  be  described  accurately  using  a  cluster  of  a 
small  number  of  atoms.  Indeed,  calculations  on  the  bulk  material  (  ) 
have  used  the  cluster  approach.  To  take  into  account  the  effect  of  the 
neighboring  bulk,  appropriate  boundary  conditions  are  imposed,  as  has 
been  discussed  in  chapter  2.  In  the  case  of  silica,  the  cluster  is 
embedded  in  a  point  charge  array,  with  the  charge  transfer  determined  in 
a  self-consistent  manner.  This  gives  the  proper  electrostatic  potential 
in  the  region  of  the  surface  cluster,  and  also  provides  for  charge 
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neutrality.  Previous  calculations  have  either  employed  a  free  space 
termination,  or  imposed  periodic  boundary  conditions  which  require 
defects  to  be  present  in  all  unit  cells  of  the  crystal. 

The  clusters  used  here  for  the  bulk  calculations  v;ere  an  SiO/j 
tetrahedron,  and  an  Si-O-Si  molecule.  Both  clusters  were  embedded  in  a 
point  charge  array  of  3*3*3  unit  cells  in  size.  The  calculations  of  the 
surface  were  done  using  a  cluster  consisting  of  the  surface  silicon 
atom,  its  two  neighboring  oxygen  atoms,  and  the  two  adjacent  silicon 
atoms  to  form  two  joined  Si-O-Si  molecules.  This  cluster  was  embedded 
in  a  point  charge  array  of  3*3*2  unit  cells.  The  position  of  the 
surface  silicon  atom  was  adjusted  in  order  to  minimize  the  total  energy. 
An  oxygen  ator.  was  then  placed  above  the  surface  silicon  in  order  to 
study  the  surface  bonding  state. 

The  one-electron  energy  levels  calculated  for  the  bulk  Si02  are  in 
good  agreement  v;ith  experiment  and  previous  theory  (fig.  ^.1).  The 
calculated  band  gap  is  9-^  eV,  as  compared  to  8.9  eV  for  experiment. 
Tne  valence  band  width  is  calculated  to  be  9.2  eV,  while  experimental 
measurements  give  about  11  eV.  The  oxygen  2s  band  is  found  at  30.3  eV 
below  the  conduction  band,  vihile  experiment  places  it  28  eV  below.  The 
valence  bands  are  found  to  be  mostly  oxygen  2p,  with  some  silicon  3s  and 
3p  character  mixed  in.  This  is  expected,  and  agrees  with  the  charge 
transfer  of  nearly  one-half  of  an  electron  from  each  silicon  to  each  of 
4  oxygen  atoms,  or  a  configuration  of  Si^+O",  as  determined  by  a 
Mulliken  population  analysis.  This  agrees  with  the  observation  that 
silica  has  both  ionic  and  covalent  properties. 
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For  the  unreconstructed  surface,  the  E'  center  is  found  to  have  a 

s 

singly  occupied  energy  level  1.7  eV  above  the  top  of  the  bulk  valence 
band  (fig.  4.2).  Stephenson  and  Binkowski  (  )  have  observed,  using 

XPS,  an  occupied  energy  level  at  0.75  eV  above  the  bulk  valence  bard 
edge  for  samples  cleaved  in  vacuum.  They  believed  this  level  to  be  an 
intrinsic  feature  of  the  bulk  Si02.  However,  their  sampling  depth  is 
only  about  30  A,  so  it  is  probable  that  this  level  is  an  E'  center  at  or 
near  the  surface  (  ).  The  method  of  sample  preparation  which  they  have 
used  (grinding)  could  easily  have  caused  these  broken  bonds  to  be  formed 
(  ).  It  is  energetically  favorable  for  an  oxygen  atom  to  bind  to  this 

surface  silicon  atom,  with  a  binding  energy  of  5.1  eV,  which  is 
considerably  greater  than  the  2.6  eV  per  oxygen  atom  necessary  to 
dissociate  molecular  oxygen.  V/e  shall  nov?  turn  our  attention  to  this 
configuration,  which  is  the  Si=0  surface  double  bond. 

The  occupied  one-electron  energy  levels  of  the  surface  state  are 

all  belov;  the  top  of  the  bulk  valence  band.  This  is  consistent  with 

optical  studies,  which  have  found  no  occupied  surface  states  in  the  band 

gap  (  ).  The  loviering  of  the  valence  band  can  be  attributed  to  the 

O-Si-0  bond  angle  at  the  surface  being  greater  than  the  perfect 
0 

tetrahedral  109  angle.  The  bottom  of  the  conduction  band  for  the 
surface  is  found  to  be  4.6  eV  below  the  bulk  conduction  band  edge.  The 
occupied  valence  surface  states  can  be  described  as  bonding  and 
non-bonding  states  between  the  adsorbed  oxygen  and  the  surface  silicon 
atom  (fig.  4.3a,b).  When  a  surface  electron  is  excited  out  of  the 
valence  band,  an  anti-bonding  orbital  (fig.  4.3c)  drops  out  of  the 
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conduction  band  and  is  singly  occupied.  Its  energy  lowers  considerably 
due  to  relaxation,  and  is  found  at  1.^  eV  above  the  top  of  the  bulk 
valence  band  (fig.  4.2).  The  position  of  the  oxygen  atom  is  0.1  a.  u. 
farther  out  from  the  surface  silicon  atom  for  this  excited  state, 
relative  to  its  ground  state  position.  This  energy  difference  is  in 
good  agreement  with  Berr.udez  and  Ritz  (  ) ,  and  the  transition  from  a 
non-bonding  to  an  anti-bonding  orbital  of  the  surface  bond  is  as  they 
have  described.  The  energy  change  associated  with  this  transition  is 
6.1  eV,  which  is  in  agreement  with  the  experimental  value  of  5.0  eV  if 
correlation  effects  are  considered.  The  7.4  eV  transition  can  be 
understood  in  terms  of  valence  band  structure,  with  transitions  from 
levels  in  the  bonding  part  of  the  valence  band  to  the  same  anti-bonding 
orbital. 

It  is  noted  that  this  model,  the  Si=0  surface  bond,  does  not 
predict  the  surface  electronic  transition  at  3.5  eV  that  is  seen  in  the 
second  derivative  ELS  spectra.  It  has  been  shown  by  Gallon  and 
Underwood  (  )  that  this  peak  is  in  fact  an  artifact  of  the  second 
derivative  mode  of  detection,  and  it  is  not  seen  in  the  non-derivative 
ELS  spectrum.  They  argue  that  the  3.5  eV  peak  is  produced  by  the 
overlap  of  the  "wings"  of  the  primary  and  the  5.0  eV  peaks.  The 
calculations  presented  here  agree  with  this  interpretation. 

One  remaining  point  of  controversy  is  the  difference  between  XPS 
and  UPS  in  discribing  the  top  of  the  valence  bands  of  Si02  (fig.  4.4) 
There  exists  at  this  time  no  agreement  as  to  the  cause  of  this 
discrepancy.  The  calculations  presented  here  do,  however,  favor  one 
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explaination.  The  two  spectra  are  spatially  sensitive  to  different 
parts  of  the  sample:  XPS  is  sensitive  to  the  bulk  for  initial  states 
with  small  binding  energy,  v;hile  40.8  eV  UPS  is  surface  sensitive.  In 
both  the  XPS  and  UPS  techniques,  the  kinetic  energy  of  the  photoelectron 
is  sufficiently  large  that  the  final  state  may  be  approximated  by  a 
plane  wave.  Therefore  the  difference  between  the  XPS  and  the  UPS 
spectra  must  reflect  the  difference  betv;een  the  bulk  and  surface  valence 
band  densities  of  states.  Comparison  between  the  XPS  spectrum  and  the 
calculated  bulk  density  of  states,  and  between  the  UPS  spectrum  and  the 
calculated  surface  density  of  states  (for  the  reconstructed  surface) 
shows  that  such  an  explanation  does  indeed  account  for  the  difference  in 
the  two  spectra  (fig.  4.4). 

The  reconstructed  silica  surface  has  been  studied  in  order  to 
determine  the  nature  of  the  low  energy  loss  ELS  peaks.  V/e  have  found 
agreement  between  the  theoretically  calculated  one-electron  energy 
levels  presented  here  and  previous  experimental  results.  Of  the  three 
models  which  can  account  for  the  peaks,  namely  a  partly  oxidized  surface 
region  SiOjj,  the  center,  and  the  Si=0  double  bond,  we  conclude  that 
the  last  is  most  likely.  Support  for  this  conclusion  comes  from 
experimental  evidence  for  reconstruction,  and  from  the  large  binding 
energy  of  the  oxygen  atom  to  the  surface. 


CONCLUSIONS 


In  this  study,  the  ab  initio  unrestricted  hartree-Fock  raethod  has 
been  used  to  calculate  the  bulk  and  surface  electronic  structure  of  of 
silica,  and  to  study  the  initial  interaction  of  oxygen  with  the  aluminum 
(100)  surface.  In  both  of  these  materials,  the  cluster  model  has  been 
shown  to  be  a  useful  tool  in  the  calculation  of  the  electronic 
structures. 

In  studying  the  interaction  af  oxygen  with  the  aluminum  (100) 
surface,  the  size  of  the  aluminum  cluster  used  had  a  considerable  effect 
on  the  results.  Since  aluminum  is  a  metal,  a  fairly  large  cluster  of 
atoms  is  needed  to  describe  the  substrate  accurately.  That  the  cluster 
model  worked  at  all  for  this  calculation  is  due  to  the  fact  that 


chemisorption  is  a  local  process 
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It  has  been  shown  that  an  oxygen  adsorbate  is  bound  most  strongly 
at  the  hole  site  of  the  clean  aluminum  surface.  Tne  equilibrium 
position  of  the  adsorbed  atom  is  about  0.2  bohr  below  the  plane  of  the 
surface-  A  charge  of  1.3  electrons  is  transfered  from  the  aluminum 
substrate  to  the  oxygen  adsorbate.  Thi'S  charge  transfer  creates  an 
electric  dipole  which  causes  the  aluminum  core  levels  to  be  shifted 
downward  in  energy  by  1.1  eV.  Previous  experiments  studying  the  work 
function,  UPS  spectrum,  and  the  EAPFS  have  reported  results  for  the 
initial  adsorption  of  oxygen  onto  this  surface  v.'hich  are  in  agreement 
with  the  results  of  the  calculations  presented  here. 

These  experimental  results  have  also  shown  that  the  nature  of  the 
adsorption  process  changes  at  about  one  nomolayer  of  coverage.  Tne 
theoretical  studies  reported  here  lead  one  to  conclude  that  at  this 
coverage  the  energetically  most  favored  sites  for  chemisorption,  the 
hole  sites,  become  fully  occupied,  and  adsorption  continues  at  the 
bridge  site.  It  is  adsorption  of  oxygen  at  this  site  which  allows  for 
incorporation  of  oxygen  beneath  the  aluminum  surface,  and  leads  to  the 
formation  of  the  oxide  AI2O2. 

Also  reported  here  is  a  study  of  the  bulk  and  surface  electronic 
structure  of  silica  (e^-quartz).  The  cluster  model  was  again  used.  The 
bulk  solid  and  the  surface  were  both  modeled  with  clusters  of  a  few 
atoms,  plus  the  appropriate  boundary  conditions.  In  this  case,  the 
boundary  condition  imposed  vjas  a  point  ion  array  to  provide  for  charge 
neutrality  and  to  provide  the  correct  Madelung  field. 
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For  the  silica  bulk,  the  calculated  one-electron  enerRy  levels  were 
compared  with  experimental  XPS  spectra,  and  the  agreement  was  found  to 
be  excellent,  as  reported  in  chapter  For  the  reconstructed  silica 
surface,  the  structure  of  the  SirO  double  bond  was  found  to  be  the 
energetically  most  favorable  surface  configuration.  Tne  electronic 
structure  of  this  surface  bond  is  able  to  explain  the  low  energy  ELS 
peaks  seen  experimentally,  lliis  calculations  reports  a  value  of  6.1  eV 
for  the  lowest  energy  transition,  as  compared  to  5.0  eV  for  the 
experimental  result. 

Finally,  the  difference  in  the  reported  valence  band  structures  as 
measured  in  XPS  and  UPS  experiments  is  explained.  Tne  XPS  method  is 
sensitive  to  the  bulk  material,  while  UPS  is  surface  sensitive. 
Comparison  of  these  spectra  with  the  calculated  valence  band  energy 
levels  shows  good  agreement  between  the  theory  presented  here  and  the 
experimenal  spectra. 

Tne  ab  initio  unrestricted  Hartree-Fock  method,  along  with  the 
cluster  model,  has  been  shown  to  accurately  describe  the  local 
electronic  properties  of  many  systems.  In  this  report,  this  method  has 
been  shown  to  be  successful  in  describing  the  adsorption  of  oxygen  onto 
the  aluminum  surface,  and  in  describing  the  bulk  and  surface  electronic 
structure  of  Si02.  Tnis  technique  has  also  been  used  to  successfully 
calculate  the  electronic  structure  of  semiconducting  polymers,  including 
defects  and  impurities  (^6),  and  in  the  description  of  the  excitonic 
structure-  of  crystalline  silicon  (^(7). 
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